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identification  of  vertical  angle  diversity  as  the  most  cost  and  perfor- 
mance effective  AAC  technique  for  tre!*ef:catter  systems. 
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development  of  a comprehensive  troposcatter  computer  model  for  pre- 
diction of  multipath  spread,  aperture  to  medium  coupling  loss,  and  angle 
diversity  parameters. 
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SECTION  1 


INTRODUCTION 

This  report  is  the  first  interim  technical  report  under  the 
Adaptive  Antenna  Control" program.  Contract  No. DAAB07-76-8-85 . The 
Adaptive  Antenna  Control  (AAC)  program  is  under  the  contract  manage- 
ment  of  the  Project  Manager  for  DCS  (Army)  Communications  Systems, 

Ft. Monmouth,  New  Jersey.  ECOM,  Ft. Monmouth,  N.J.  is  responsible  for 
technical  direction  of  the  program.  SIGNATRCN,  Inc.  of  Lexington,  MA 
with  RF  Systems,  Inc.  of  Cohassett  as  a major  siibcontractor  have 
the  responsibility  for  executing  the  program  tasks.  This  report 
describes  the  AAC  program,  summarizes  system  analysis  work  to  date, 
and  projects  the  scope  and  direction  of  future  efforts.  Hardware 
descriptions  and  specific  design  information  for  the  implementa- 
tion tasks  will  be  provided  in  the  AAC  design  report  due  22  January 
1977. 


1.1  Need  for  Adaptive  Antenna  Control 

Transhorizon  radio  systems  within  the  Defense  Communications 
System  (DCS)  employ  both  tropospheric  scatter  and  diffraction  over 
obstacles  as  primary  propagation  modes.  Troposcatter  systems  suf- 
fer from  short  term  multipath  fading  and  long  term  depressions  in 
median  received  signal  level.  The  short  term  fading  can  be  al- 
leviated by  the  use  of  additional  diversity.  Angle  of  arrival 
diversity  represents  an  attractive  method  of  realizing  this  addi- 
tional diversity  as  it  requires  no  additional  frequencies  or 
antennas.  Depression  of  the  median  received  signal  is  primarily 
due  to  substandard  refractive  index  gradients  which  bend  the  radio 
wave  upward  and  increase  the  scattering  angle  and  hence  the  path 
loss.  This  condition  cannot  be  alleviated  by  a change  in  antenna 
pointing  angle  as  the  change  in  refractive  index  causes  more  up- 
ward bending  of  the  radio  wave  but  maintains  the  same  take-off 
angle.  Note  that  downward  bending  of  the  radio  wave  resulting 
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from  superrefractive  or  ducting  conditions  is  not  a problem  on 
tropcscatter  circuits  because  any  decrease  in  signal  due  to  earth 
blockage  is  more  than  compensated  for  by  the  decreased  path  loss 
due  to  a decrease  in  the  scatter  angle. 

For  transhorizon  radio  links  using  diffraction  over  an  obstacle 
(3cnife-edge)  , superrefractive  conditions  are  the  primary  cause  of 
signal  outages.  A knife-edge  diffraction  path  is  characterized 
by  either  the  absenc;e  of  short  term  fading  due  to  scattering  or 
milder  forms  of  short  term  fading  due  to  multiple  path  interference 
conditions.  Thus  outages  tend  to  be  long  term  and  are  primarily 
produced  by  downward  bending  of  the  radio  wave  and  resulting  earth 
blockage  during  superrefractive  conditions. 

Experimental  tests  on  a DCS  diffraction  system  have  shown  these 
outages  to  exist  and  because  of  their  long  term  nature,  they  rep- 
resent a serious  deterioration  in  link  availability  . Clearly 
an  automatic  method  of  varying  the  transmit  and  receive  antenna 
beams  in  elevation  could  provide  protection  against  these  outage 
conditions.  On  the  other  hand  the  absence  of  a significant  multi- 
path  fading  phenomenon  on  most  diffraction  links  suggest  that  addi- 
tional diversity  in  the  form  of  angle  diversity  is  not  of  appreciable 
value  provided  adaptive  antenna  pointing  can  be  accomplished. 

The  DCS  tropospheric  scatter  radio  systems  may  be  faced  with 
requirements  to  operate  at  higher  carrier  frequencies  in  certain 
host  countries.  Normally,  for  wide  beam  antennas  and  short  paths, 
the  additional  antenna  gain  at  a higher  frequency  more  than  off- 
sets the  increased  path  transmission  loss.  However,  for  longer 
paths  and  generally  narrower  bearawidths,  at  higher  frequencies 
the  aperture  to  medium  coupling  loss  will  reduce  the  effective 
antenna  gain  and  a net  loss  will  result  after  frequency  conversion. 
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Adaptive  antenna  techniques  offer  the  potential  of  mitigating 
coupling  loss  and  meeting  system  performance  requirements  at 
higher  frequencies  on  existing  links. 

1 • 2 AAC  Approach  for  DCS  Application 

As  a result  of  the  above  considerations,  an  adaptive  antenna 
control  system  has  been  proposed  which  utilizes  angle  diversity  on 
DCS  troposcatter  systems  and  antenna  beam  steering  on  DCS  diffrac- 
tion systems.  The  angle  diversity  system  will  employ  dual 
polarized  feedhorns  to  produce  vertically  displaced  antenna 
beams  which  intersect  at  or  near  the  half  power  beamv;idth  points. 

In  this  manner,  the  scattering  angle  and  hence  the  path  loss 
of  the  elevated  beam  is  minimized  while  maintaining  a significant 
decorrelation  between  the  received  diversity  signals.  For  example  a 
168  mile  path  with  28 ' antennas  at  4 . 46  GHz  has  a relative  loss  of  1 . 7 dD  in 
the  elevated  beam  and  power  correlation  coefficient  of  0,52.  An 
analysis  shows  that  almost  all  DCS  troposcatter  systems  will 
realize  significant  improvement  from  the  addition  of  angle  diver- 
sity. In  some  cases  the  additional  angle  of  arrival  diversity  may 
complement  the  existing  quadruple  diversity  operation  to  produce 
approximately  8th  order  diversity  in  a dual  space/dual  frequency/ 
dual  angle  configuration.  In  other  applications  dual  frequency 
diversity  systems  may  be  upgraded  to  quad  diversity  perfor- 
mance by  using  a dual  frequency/dual  angle  configuration.  A re- 
duction to  this  configuration  from  a conventional  duplex  four- 
frequency,  four-antenna  quad  diversity  results  in  the  savings 
of  two  antennas  which  might  prove  to  be  an  attractive  alternative 
for  certain  links.  Similarly  a reduction  to  dual  space/dual  angle 
provides  nearly  quadruple  diversity  performance  and  saves  two  fre- 
quencies . 

Performance  predictions  for  classes  of  DCS  links  must  be  performed 
in  the  study  as  angle  diversity  (AD)  does  not  yield  an  exact  doubling  of  ^ 
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the  diversity  order  due  to  loss  in  the  elevated  beam  and  correlation 
between  diversity  signals.  The  resulting  performance  improvement 
for  defined  classes  and  the  cost  analysis  for  AD  implementation  will 
determine  the  cost  effectiveness  of  this  approach  for  each  class. 

For  DCS  diffraction  systems,  a mechanical  means  of  providing 
antenna  movement  and  an  electronic  adaptive  control  system  for 
optimally  pointing  the  antenna  is  proposed.  Electronic  steering 
has  been  judged  cost  ineffective  because  of  the  difficulty  of 
providing  amplitude  and/or  phase  control  of  large  transmitter 
powers  (~10  kw)  . A meciianiaal  design  of  a weatherized  moveable 
feedhorn  structure  is  under  development  by  RF  Systems,  The  AAC 
program  will  include  the  design  of  an  adaptive  antenna  aiming  sys- 
tem and  an  analysis  thereof. 

Within  the  AAC  program,  which  is  developing  an  angle  diversity 
feed  system  and  the  application  guidelines  for  using  angle  diversity 
to  enhance  the  performance  of  digital  troposcatter  in  a frequency 
congested  environment,  it  is  desirable  and  convenient  to  also  address 
the  effects  of  aperture  meditim  of  coupling  loss  in  the  S and  C 
frequency  ban<i  in  order  to  extend  the  usable  range  of  troposcatter 
systems  operating  in  these  bands  (e.g.,  TRC-170) . Under  the  subject 
program,  troposcatter  tests  to  be  performed  at  RADC  for  angle  diver- 
sity present  an  opportunity  to  obtain  empirical  data  to  extend  the 
theoretical  model  of  coupling  loss  to  higher  frequencies  and  de- 
velop candidate  techniques  capable  of  mitigating  coupling  loss.  Be- 
cause of  the  likelihood  that  angle  diversity  will  be  used  in  the  DCS, 
the  suitability  of  angle  diversity  versus  other  candidate  techniques 
such  as  multiple  apertures  also  will  be  addressed  to  determine 
whether  angle  diversity  is  superior  to  other  techniques  in  a joint 
cost  performance  tradeoff.  As  such,  calculations  for  DCS  links  are 
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required  to  determine  the  resultant  change  in  required  system 
gain  and  system  cost  after  frequency  conversion  and  the  introduc- 
tion of  each  candidate  technique  in  order  identify  the  optimum 
improvement  technique. 

1.3  AAC  Advanced  Development  Model 

The  Advanced  Development  Model  (ADM)  includes  the  design  of 
an  adaptive  antenna  control  system  with  mechanical  movement  for  dif- 
fraction systems  and  the  implementation  of  an  angle  diversity  sub- 
system for  troposcatter  systems.  The  adaptive  antenna  aiming  sys- 
tem utilizes  a servomechanism  for  antenna  control  and  a straight- 
forward adaptation  algorithm  to  maximize  the  received  power  at  each 
link  terminal.  It  is  envisioned  that  control  information  would  be 
transmitted  over  an  order  wire  channel  primarily  to  preclude  both 
antannas  from  trying  to  adapt  simultaneously. 

The  ADM  for  a troposcatter  system  consists  of  a vertical  ar- 
ray of  two  dual  polarized  feedhorns,  additional  waveguide  and  RF 
down  converters  for  each  additional  diversity,  and  a predetection 
combiner  for  digital  channel  applications.  The  PDC  combines  up  to 
four  digital  angle  diversity  signals  and  considers  the  presence 
of  multipath  by  weighting  each  angle  diversity  signal  with  a three 
tap  transversal  filter.  The  output  of  the  PDC  is  combined  with 
the  digital  modem  combination  of  the  dual  space/dual  frequency 
diversity  channels  to  provide  up  to  eight  order  diversity.  The 
PDC  will  be  designed  for  possible  operation  with  any  linear  modula- 
tion modem  type,  e.g.,  2 PSK,  QPSK,  Staggered  QPSK,MSK,  and  provides 
optimum  combining  when  the  modem  interfaces  the  PDC  with  a one  bit 
error  signal.  Although  this  error  signal  interface  should  be  pos- 
sible for  any  of  the  linear  modulation  modem  types,  the  PDC  can  be 
configured  in  a suboptimum  mode  for  combining  without  any  inter- 
face signals  between  the  modem  and  PDC. 
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The  PDC  concept  is  based  on  technology  developed  during  the 
Megabit  Digital  Troposcatter  Subsystem  (MDTS)  program  and  the  PDC 
design  includes  4 card  types  identical  to  those  used  in  the  MDTS 
modem.  Interface  and  excellent  operation  of  the  PDC  with  the  MDTS 
modem  is  guaranteed  by  the  functional  and  design  similarity  betv/een 
the  two  systems. 

1.4  AAC  Field  Test 

The  PADC  test  bed  (Youngstown-Verona,NY)  to  be  used  in  the  field 
test  is  a 168  mile  tropospheric  scatter  radio  system.  We  propose 
to  fabricate  a simplex  angle  diversity  system  using  tv/o  vertical 
feedhorns  at  each  of  the  two  receiver  (Verona)  test  link  antennas. 
The  maximum  diversity  complement  achievable  under  this  configura- 
tion will  be  eighth  order  for  a dual  space/dual  frequency/dual 
angle  (2S/2P/2A)  system.  A Predetection  Combiner  will  be  fabricated 
to  provide  the  diversity  combining  function  for  the  digital  tests. 

It  is  planned  that  three  MDTS  modems  will  be  supplied  for  the 
digital  tests.  For  a 2S/2P/2.A  configuration,  one  modem  is  required 
for  transmission*,  one  modem  is  required  to  combine  the  conventional 
quad  diversity  channels  while  the  PDC  combines  the  four  additional 
angle  diversity  channels,  and  the  third  modem  is  used  at  the  re- 
ceiver for  a simultaneous  comparison  test  of  quad  diversity  with- 
out angle  diversity. 

In  addition  to  the  digital  transmission  tests,  extensive 
testing  of  channel  characteristics  such  as  multipath  spread,  diver- 
sity channel  correlation,  aperture-to-medium  coupling  loss,  and  re- 
ceiver power  statistics  is  planned.  These  tests  will  serve  to  vali- 
date the  performance  model  of  the  AAC  Advanced  Development  Model 
and  they  will  also  be  used  to  improve  and  validate  the  troposcatter 


* This  modem  can  be  an  MDTS  modulator  since  reception  is  not 
required  at  Youngstown, 
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channel  model.  Of  particular  interest  is  an  accuract  prediction 
model  for  multipath  spread  and  coupling  loss  which  includes  an 
assessment  of  the  variability  of  these  quantities  as  well  as  the 
median  values.  An  examination  has  also  begun  to  determine  the 
^^^^i^^ility  of  meteorological  data  which  could  be  used  to  relate 
atmospheric  measurements  to  communication  parameter  prediction. 

Some  of  the  primary  goals  of  the  test  program  include 

• verification  of  the  channel  characteristics  predic- 
tion model 

» verification  of  the  angle  diversity  performance 
prediction  model 

• measurement  of  channel  parameters  including  aperature- 
to-medium  coupling  loss 

• establishing  the  operational  capability  of  an  AD 
troposcatter  system 

» determination  of  equipment  parameter  modifications 
in  the  AAC  model  to  improve  performance  and  cost 
effectiveness 

• establishing  the  interface  capability  of  the  AAC 
model  with  digital  equipments. 

The  schedule  for  the  accomplishment  of  the  AAC  tasks  and  field 
tests  outlined  in  this  section  is  shown  in  Fig.  1.1. 

1 . 5 Report  Summary 

A theoretical  model  for  troposcatter  communication  has  been 
developed  as  the  basis  for  prediction  of  multipath  spread,  angle 
diversity  communication  parameters,  and  aperture— to— medium  coupling 
loss.  The  model  derives  from  Tartarskii's  original  work  on  scat- 
tered signals  due  to  random  variations  in  the  refractive  index  in 
the  common  volume  defined  by  the  antenna  beams.  Under  a decorrela- 
tion assximption  for  the  refractive  index  correlation  function  and  a 
limitation  on  the  signal  bandwidth,  the  uncorrelated  scattering 
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Gaussian  conmiunication  model  can  be  used  to  represent  signal  transfer 
over  troposcatter  paths.  The  result  is  that  a set  of  L diversity  chan- 
nels is  completely  characterized  by  the  set  of  cross-channel  multi- 
path  profiles  defined  by 

h (t)h  (t)  = 6(t-T)  Q.„(t) 
t m tm 


where  and  are  the  random  troposcatter  impulse  responses, 

for  the  /,th  and  mth  diversity#  respectively.  For  an  angle 
diversity  system  where  there  exists  significant  correlation  only 
between  the  main  and  squinted  beam#  three  profiles;  for  the 

main  beam,  for  the  squinted  beam,  and  the  product  beam, 

are  required.  Section  2 presents  the  theoretical  basis  and  the 
volxome  integration  formulation  for  the  determination  of  the  profile 
function  Q^j^j(t)  . Also  in  this  section,  formulas  and  asymptotic 
results  for  multipath  spread  and  coupling  loss  are  derived. 


Accurate  prediction  of  channel  parameters  in  accordance  with 
the  theory  of  Section  2 requires  an  integration  over  the  scattering 
volume  of  the  antenna  patterns  and  the  refractive  index  spectrum. 

An  exact  calculation  method  for  a specific  refraction  profile  is 
used  as  the  basis  for  the  numerical  integration  procedure.  A de- 
scription of  the  refraction  profile  model  and  the  integration 
geometry  begins  Section  3.  Following  this  description  we  present 
and  interpret  computation  results  for  multipath  spread  and  coupling 
loss.  Much  of  the  discrepancy  between  predicted  and  measured  multi- 
path  spread  is  shown  to  a result  of  measurement  inaccuracy.  Trans- 
mitter filter  effects  and  prober  processing  inaccuracies  are  the 
most  significant  contributors  and  their  effects  can  be  dramatic  when 
measuring  small  multipath  widths  (2a  < 100  nanoseconds) . The  MDTS 
field  test  data  was  used  to  compare  predicted  and  measured  multipath 
spreads  and  the  results  of  Section  3 serve  to  validate  the  utility 
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of  the  prediction  model  used.  Following  the  multipath  spread 
calculations  and  data  analysis,  a similar  set  of  coupling  loss 
computations  and  comparison  with  previous  empirical  results  has 
been  initiated  and  described  in  Section  3.  For  zero  take-off 
angles  of  the  antenna  beams  the  volume  integration  method  generally 
predicts  larger  coupling  loss  then  available  semi-empirical  tech- 
niques. Better  agreement  results  when  the  volume  integration 
method  assumes  positive  take-off  angles.  Antenna  take-off  angles 
(the  angle  between  the  antenna  boresight  and  the  local  horizon) 
are  shown  to  be  a significant  link  parameter  which  needs  to  be  ac- 
curately specified  for  useful  performance  prediction.  In  the 
absence  of  this  information  it  is  not  clear  whether  the  prediction 
model  is  imprecise  or  incorrect  link  parameter  specification  is 
responsible  for  the  inaccuracies.  Also  in  this  section  an  assess- 
ment of  multi-feedhorns  (angle  diversity)  vs.  multi-antennas  (space 
diversity)  for  the  mitigation  of  aperture-to-medium  coupling  loss 
is  presented.  The  angle  diversity  approach  is  shown  to  have  a 
clear  superiority  from  a cost/performance  standpoint. 

In  Section  4 the  methods  of  adaptive  antenna  control  are  in- 
troduced and  evaluated.  An  adaptive  mechanical  feedhorn  arrangement 
is  recommended  for  diffraction  links  and  a dual  vertical  feedhorn 
with  angle  diversity  combining  is  recommended  for  troposcatter  links 
An  evaluation  of  the  potential  utility  of  angle  diversity  in  the  DCS 
is  provided.  As  an  approximate  guide,  angle  diversity  produces  a 
significant  improvement  if  the  ratio  of  beamwidth  to  scattering  angl 
is  less  than  0.875.  This  criterion  is  satisfied  by  most  DCS  tropo- 
scatter systems.  This  section  concludes  with  a summary  of  previous 
adaptive  antenna  control  experiments. 
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Specific  implementation  factors  for  adaptive  antenna  control 
on  both  diffraction  and  troposcatter  links  are  presented  in  Sec- 
tion 5.  Diversity  combiner  structures  for  both  analog  and  digital 
applications  are  also  discussed.  In  Section  6,  performance  improve- 
ment considerations  are  addressed  for  the  diffraction  path.  Varia- 
tions in  take-off  angle  due  to  fluctuations  in  the  refractive  index 
gradient  are  proportional  to  path  distance  and  independent  of  ob- 
stacle height.  For  a 175  nautical  mile  path,  realistic  refractive 
index  variations  could  produce  a change  of  + ^5°.  Thus  a fixed 
antenna  beam  system  with  a beamwidth  less  than  or  equal  to  ^5°  could 
experience  a complete  loss  of  signal.  The  performance  of  angle 
diversity  with  analog  transmission  techniques  is  computed  for  the 
RADC  test  link.  Finally,  digital  performance  factors  are  considered 
and  the  proposed  method  of  computing  modem  performance  is  indicated. 

A major  portion  of  the  study  program  within  the  next  reporting  period 
will  entail  the  development  of  modem  performance  characteristics  for 
typical  and  projected  angle  diversity  links. 
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SECTION  2 

TROPOSCATTER  COMMUNICATION  MODEL 

In  this  section  we  develop  the  theoretical  basis  for  a tropo- 
scatter  coiranunication  model  which  will  be  used  to  relate  atmospheric 
parameters  to  communication  parameters  and  to  relate  system  perfor- 
mance to  the  communications  parameters.  The  atmospheric  parameters 
of  interest  are  those  available  in  conventional  radio-sonde  mea- 
surement techniques,  i.e.,  air  pressure,  temperature,  and  water  vapor 
pressure  as  a function  of  height.  These  parameters  can  in  turn  be 
related  to  the  statistical  description  of  the  refractive  index 
fluctuations.  This  statistical  description  is  used  to  develop  the 
communication  parameters  of  an  equivalent  random  linear  filter. 

This  section  concerns  itself  with  the  statistical  description  of 
this  filter.  Numerical  computation  of  communication  parameters 
such  as  multipa'-h  spread  and  antenna  coupling  loss  are  presented 
in  Section  3 . In  Section  6 the  methods  of  using  these  parameters 
to  compute  communication  performance  is  addressed.  The  relationship 
between  atmospheric  and  communication  parameters  and  the  techniques 
for  obtaining  meteorological  information  will  be  dealt  with  in 
greater  detail  in  subsequent  reports.  Although  the  communications 
model  is  applicable  to  both  analog  and  digital  systems,  emphasis  is 
placed  on  the  digital  system  application  since  the  multipath  effects 
are  more  pronounced  on  the  generally  wider  bandwidth  digital  systems . 

2 . 1 Uncorrelated  Scatter  (US)  Channel 

The  troposcatter  channel,  viewed  as  a black  box,  has  an  input- 
output  relationship  identical  to  a randomly  time-varying  linear 
filter.  Because  the  time  variation  is  very  slow  (~  1 Hz)  compared 
to  the  large  signal  bandwidths  (1-20  MHz) , performance  can  be  deter- 
mined from  an  ensemble  description  of  an  equivalent  random  linear 
filter.  An  ensemble  member  characterization  is  shown  in  Fig.  2.1 


2-1 


in  terms  of  fixed  tran;5mit  and  receive  filters  and  a random  filter 
representing  the  troposcatter  medivim. 

* 

The  composite  filter  is  the  threefold  convolution 


h^(t)  = h^(t)  0f^(t)  ©hj^(t)  (2.1) 

where  the  symibol  © denotes  convolution.  If  we  define  the  coiribined 
fixed  filter 


then 


= h^(t)  * hj^(t) 


I 

h (t-u)  f , (u)  du. 

XK  V 


The  impulse  response  f^ (t)  represents  the  random  tropospheric 
scatter  channel.  An  important  communications  model  for  this  chan- 
nel is  the  complex  Gaussian  uncorrelated  scatter  (US)  channel  [2.1]. 
This  model  defines  zero-mean  Gaussian  with  uncor- 

related scattering  for  different  transmission  delays,  i.e., 

f^(t)  fjjj*(f)  = 6(t-t’)  (2.2) 

where  defined  as  the  cross  channel  multipath  profile.  Note 

that  for  this  model  the  composite  filter  has  second  moment 


■ ' — - -■  OD  00 

‘ I J ‘’to 


' h^(t-u)  h^(f-u)  Q^^(u)  du. 


Complex  notation  to  represent  real  bandpass  signals  is  used  through- 
out the  report.  The  real  impulse  response  is  Re{h.(,(t)  where 

±0  is  the  carrier  frequency  which  is  assumed  to  be  much  larger  than 
the  signal  bandwidth. 
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If  the  fixed  filters  have  a maximvun  bandwidth  of  B Hz,  a 
sampled  data  representation  of  the  composite  filter  is 


h^  (t)  = Z h^^  sine  (t/i-i) 


where 


T < 1/B 


li^  = ^ sine  (t/T-i)  dt  = 


and 


sine  X = sin(TTx)/^x. 


From  Nyquist's  Theorem,  the  samples  h^^  are  a complete  rep- 
resentation of  the  multichannel  filter.  A complete  statistical 
description  of  the  communication  channel  is  thus  obtained  if  h^^ 
is  described  statistically.  A complete  description  for  the  special 
case  of  zero  mean  Gaussian  samples  is  given  by  the  second  moment 

h.,h.  . For  the  US  model,  this  second  moment  is 
1-0  3m 

' 00 

and  for  typical  transmitter  and  receiver  filters  which  do  not  sig- 
nificantly distort  the  information  signal,  the  width  of  h (*)is 
such  that  when  i ^ j the  value  of  the  correlation  function  is  much 
smaller  than  when  i = j . Thus  to  a good  approximation 


h.,h.  = a.r  lh^(iT-u)l  Q,  (u)  du. 

i-t  ^m  13J  TR  ' -tm 


We  will  now  show  that  for  Gaussian  distributed  refractive  index 
fluctuations,  if  a condition  on  the  decorrelation  of  these  fluctua- 
tions is  met,  the  troposcatter  channel  is  well  approximated  by  the 

* A two  sided  bandwidth  B is  used  with  complex  notation  corresponding 
to  the  RF  bandv/idth. 
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uncorrelated  scatter  channel  and 

profile  Q,  (t)  is  expressible  as 
-Cm 

to  the  combined  antenna  patterns 
index  f lv.ctuations . 


the  cross  channel  multipath 
a volvime  integral  with  respect 
and  the  spectrum  of  the  refractive 


Consider  transmission  through  a particular  troposcatter  chan- 
nel with  the  geometry  shown  in  Fig.  2.2.  At  a scattering  point 
r ' , which  can  be  considered  to  be  in  the  far  field  with  respect 
to  both  the  transmitter  (T)  and  receiver  (R) , we  assume  an  incident 

electric  field  A (r ' ) — which  is  a plane  wave  with  orogration  vector 

o ~ ■ ■ 

>:  m where  k = 2Trf^/c  =»  2rr/X , f^  = carrier  frequency,  c = speed  of 
light  in  the  medium,  and  X = wavelength.  A fraction  of  the  energy 
associated  with  the  incident  field  is  scattered  in  the  direction  n 
toward  the  receiver.  The  scattering  mechanism  can  be  attributed 
to  random  spatial  variations  in  the  refractive  index  at  the  point 
r ' . Letn^(r')  be  the  random  component  of  the  refractive  index, viz., 

n(r' ) = n(r' ) + n^  (r' ) , 


then  in  the  Appendix  following  we  show  a derivation  first  developed  by 
Tatarski  [2.2.]  that  the  electric  field  at  the  receiver  antenna  is 
given  by 


e(r) 


jk( Ir-r 

e-'  


+ m*r ' ) 


(2.4) 


Troposcatter  antennas  utilize  a receiving  aperture  typically 
in  the  form  of  a parabolic  reflector  with  a feedhcrn  placed  at  the 
focal  point.  The  mean  field  at  the  receiver  is  obtained  by  in- 
tegrating over  the  receiving  aperture 


E(r)  = 2 JJ  e(r  + u)  d§dT|» 
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] 

i 


where  u is  in  the  plane  of  the  aperture  and  points  to  the  dif- 
ferential surface  element  dPdT].  Since  r » u,  the  effect  of 
changing  r solely  impacts  the  exponent  in  (2.4).  A useful  Taylor 
series  expansion  is 


II  1 1 r ^ 

X+U  = X + — x*u+~  u 

I _ _ 2x  L 


(u  • u) 


21  + 


and  if 


1 2 2 
kU  _ TTU 

2x  Xx 


« n. 


the  exponential  term  is  well  approximated  by 

Ml  ^ I j^(x  + — X • u) 
ik'x  + u . X — — 

e-<  _ g 


(2.5) 


that 


If  D is  the  receiving  aperture  diameter,  it  is  easy  to  verify 


« 1 


for  practical  troposcatter  pareuneters.  The  mean  field  is  thus 


. k P,  / , V , jT^(|r  - r'|+  m * r'),3  , - , . 

E(r)  = ’ )n^(r’ ) e-^  i - - dr'-f^(n) 


where 


f,(n)  = ^ II  e-^lcn  -u 


(2.6) 


(2.7) 


is  the  normalized  antenna  gain  pattern  at  the  receiver.  In  (2.6) 
we  have  suppressed  the  time  dependence  in  the  notation.  Let  u(x'»|) 
represent  a three  dimensional  delta  function  for  the  delay  £ from  the 
transmit  antenna  to  the  source  point  r'  and  down  to  the  receive  antenna,  I 


i.e. , 


U(r',  f)  = 6[U’  - sj  + |r  - r'l  - c?]. 
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The  received  field  at  time  t,  E(r,  t) , for  an  excitation  A (r',t)e^-*-' 

o — 

at  the  source  point  is  the  superposition  integral 

,2  eo 

= -^  J J U(r',5)A^(r',t-^)f^(n)n^(r')e^’^(lr-r'l+m-r')  3 , . 

V 

The  time  dependence  of  n^(r')  is  ignored  since  the  refractive 
index  does  not  change  over  the  range  of  differential  delays  of  sig- 
nificance in  the  inside  integral.  Since  the  system  is  linear  we  can 
combine  the  transmitter  and  receiver  filters  in  the  expression  for 


A . Thus  we  define 
o 

A (r',t)  = 
o 


where 


T^T 


4ns 


2 


(2.8) 


T 

s 

fy(Sl) 


transmitted  power 

transmitter  antenna  maximum  power  gain 

mean  distance  from  transmitter  to  common  volume 
normalized  transmitter  antenna  pattern  (see  Eq.(2.7)) 

h (t)  (*)h  (t)  = combined  fixed  filter  response. 

T R 


In  an  angle  diversity  system  with  one  transmit  beam  and  multiple 
receive  beams,  the  mean  field  for  the  'tth  diversity  is  given  by 


E.  (r,t)=:|^  A-?  f h„„(t-?)d? 


2Trr  I 2 J TR 
4tts  -® 


O j. 

V 

We  have  ignored  the  slight  difference  between  r and  the  actual 
distance  to  an  elevated  common  volume  in  the  denominator  of  (2.9) . 
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Also,  the  path  length  d is  approximately 


d = r + s.  (2.10) 

Since  (2,9)  represents  the  mean  field  at  the  receiver  for  im- 
pulse function  excitation,  we  must  multiply  by  the  antenna  aperature 
area  to  obtain  the  sampled  data  impulse  response.  Since  the  apera- 
ture area  is  related  to  the  receive  antenna  maximum  power  gain  by 
2 

= 4tt£A  / we  have 

where  G is  the  receive  antenna  maximum  power  gain.  For  convenience 
R 

let 

P„G  G 
T T R 

^o  ~ — 4tt — ~ free  space  received  power  per  unit  distance 

g, ( r')  = f (m)f, (n)  = normalized  combined  antenna  pattern 
■t/  “ o ~ “ 

and  define  the  smoothed  volume  integration  in  (2.9)  by 

CO 

1 1 “<£'•?)  = r 

mm  CO 

V V. 

1 


with  these  definitions,  the  sampled  data  impulse  response  is 


^ J^o  r , . , , . * ik  ( I r-r ' I +m  • r ' ) ,3  , 

_ — — I g (r* ) n (r ' ) e"^  ' ' — — dr'. 

2Vnrs  ^ 

i 


(2.13) 


The  smoothed  integration  defined  in  (2.12)  can  be  viewed  as 
the  weighted  svim  of  scatter  points  which  have  transmission  delays 
which  roughly  fall  in  the  interval  [iT-r/2, iT+x/2] . The  weighting 
falls  off  as  the  transmission  delay  moves  away  from  the  midpoint  ir . 

The  refractive  index  fluctuations  can  be  reasonably  approxi- 
mated by  a stationary  zero-mean  complex  Gaussian  process.  The 
statistics  of  h. . are  then  completely  specified  in  terms  of  the 
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refractive  index  correlation  function 


n^(r'+^)  n^*(r')  = «?  {^)  (2.14) 

or  its  three  dimensional  Fourier  transform,  the  refractive  index 
spectrum 


§ (k)  = Icp  (£_)  e“^—  ^ d^p  . 
8rr 


(2.15) 


V 


The  correlation  radius  is  defined  as  the  radius  beyond  which 
cp(£)  is  small  compared  to  its  maximxim  value. 


h . , h . = 

xl  jm 


where 


The  impulse  response  correlation  function  is 

JJ  g^(^.')gjj^(£“)cp(^"“^' d^r'd^r"  (2.16) 


k^P 


2 2 

4tt  r s 


V.V. 
1 1 


O'®  ~ “ £."1  + I5i(t')  •£'  - ni(r")  *r " . 


In  order  to  reduce  (2.16)  further,  we  must  either  know  the 
correlation  function  cp(r"-r')  or  some  estimate  of  its  correlation 
radius.  We  will  show  that  the  basis  for  the  widely  used  uncor- 
related scattering  assumption  is  that  the  correlation  radius  sat- 
isfies 

t <<  c/B8  uncorrelated  scatter  assumption  (2.17) 

o 

where  is  the  scattering  angle  at  the  lowest  point  in  the  common 
volume . 


Since  S « 1 it  follows  that  on  the  order  of  <c/B  satisfies 
the  constraint,  e.g.,  for  a bandwidth  of  10  MHz,  the  correlation 
radius  must  be  on  the  order  of  100  feet  or  less. 
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The  smoothed  integrals 


1 J = J J “(4  ’l>  I \r<:'^-?2>"52  I 


V.  V.  -» 
1 1 


V 


V 


reduce  to  two  volume  integrals 


I 1 “ v("  - ='(£'>/=)  J V I'jf  - ' 

V.  V.  V 

1 D 

where  x is  the  path  distance  through  the  source  point/ 
x(r')  = j_r'  -^1+  1 • 

Introducing  the  change  of  variables 
p = r"  - r' 

we  have 


I I ' 1 ■ ="'''>/=)  I \r  - ='<£'  p'/=) 


V.  V.  V 

1 j 


V 


The  distance  parameter 


x(r'+_p)=  jr'  + £ - s|  + If  - _P  - f"l 

can  be  reduced  by  using  the  expansion 

~~2  2 
+2£*u  + u 

2 r • u 

r / 1 + : — * r + a • u 

^ 2 — r — 

r 

= r/r  . 


where 


a 

— r 
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u (r"  / 


One  finds 

x(r'  + p)  - + m(r')  • p + jr  - r' j - n(r')  • p 

= X (r ' ) + (in  - n)  * p 

The  second  term  is  very  small  compared  to  the  first  and  is  of 
the  order 

jm-nj  1 p|  = p sin  G(r ' ) * p 6 (r ' ) * p 6 # 

where  Q(r')  and  9 are  the  scattering  angles  at  the  source  point  r' 
and  at  the  lo\7est  point  of  the  common  volume,  respectively. 

Returning  to  the  impulse  correlation  function,  we  have 
* 

J - *(£')/=X(£')a  r' 

■ J ’'TR(3f-’=(£:')/=-8p/=)  9^*  (r'  + p)'!>(p>e^^'‘ 

V - - - d^p 

Since  the  correlation  distance  is  on  the  order  of  feet  and  the 
antenna  beam  pattern  spans  distances  on  the  order  of  miles,  it 
follows  that  in  the  second  integral  where  ] p|<  -(/q  can  mahe  the 
approximation 

Also  by  using  a Taylor  series  expansion  as  in  (2.4) , and  since 

2 

/\r  « 1,  we  find  we  can  reduce  the  exponent  a as  follows.  Since 
Ir  - r'  - ^1  = 1^  - r'l  ” « * _P 

and 

m(r'  + £)  • (r'  + £)  = m * (r'  + £)  , 
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we  have 


a = m(3:')  • r'  - m(r'  + p)  • (r'  + p) 


a ~ “H  * p in  • p 


The  correlation  function  bejcomes 


h.  h.  = 

i-t  ini  .2 

4Tir  s 


■J  (p' (2.18) 

e^p 

where  x,  m,  and  n are  functions  of  the  source  point  r' . Note  in 
(2.18)  that  the  second  integral  is  a smoothed  3-dimensional  Fourier 
transform  similar  to  the  refractive  index  spectrum  defined  in  (2.15) 

Now  the  uncorrelated  scattering  assumption  states  that  the 
correlation  radius  satisfies 

« c/B8  , 

hence  in  the  second  integral  of  (2.18)  the  terminal  equipment 
impulse  response  varies  over  the  delay  range 

1 

c B 

where  B is  the  bandwidth.  Since  h (•)  is  approximately  constant 

xR 

over  an  interval  small  compared  to  a reciprocal  bandwidth  it  can 
be  taken  outside  of  the  second  integral  with  the  result 
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h.  h. 
li.  jm 


2 2 

- V 


, § k(m  - n)^  d r'  . (2.19) 


If  we  define 


Q (t)  = j U(r',  t)g  (r')g  (r’ ) $rk(m  - n)  I d‘'r 

rm  ^~m~L  J 


,1  ..3^, 


we  can  express  the  correlation  function  as 

2 2 

VV*°”2^  I 

r s -® 


However,  since  (2.21)  is  identical  in  form  to  (2.3),  it 

follows  that  0 (t)  of  (2.20)  is  the  cross  channel  multipath  pro- 
^in 

file  defined  in  (2.2)  for  the  ancorrelated  scatter  channel. 


Thus  we  have  established  that  for  a refractive  index  correla- 
tion radius  less  than  the  ratio  of  the  speed  of  light  to  the  signal 
bandwidth,  the  uncorrelated  scatter  model  is  a good  representation 
for  troposcatter  communication.  Under  ducting  conditions  or  the 
presence  of  large  atmospheric  layers,  there  will  be  more  correla- 
tion between  adjacent  delay  cells  with  the  result  that  intersym- 
bol interference  will  not  be  as  severe  and  implicit  diversity 
gains  will  not  be  as  large.  However  the  ducting  and  large  layer 
are  typically  associated  with  stronger  received  signals  and  thus 
are  of  less  importance  in  predicting  communication  outage  per- 
formance. Thus  for  the  performance  predictions  to  be  developed 
in  this  effort,  the  US  model  with  a cross  channel  multipath  pro- 
file defined  by  (2.20)  will  be  used  throughout. 
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2.2  Troposcatter  Parameters 


1 


The  complex  Gaussian  uncorrelated  scattering  model  for  tropo' 
scatter  communication  is  completely  defined  by  the  cross-channel 
multipath  profile 


f (t)f  *(t')  = 6(t-f)  Q,^(t) 

4 <•  -f,m 

which  was  determined  in  the  previous  sub-section  to  have  the  fom 


(t) 


2 2 

2tt  k P 


2 2 
r s 


f 

t 

J 

V 


•U(r',  t)  g (r‘)  g*  (r')  § Ck(m-n)]  d^f*  (2.22) 
-c,  m — 


where 

k = 2n/X,  X = wavelength 

Pq  = free  space  received  power  per  unit  distance 

r s=  mean  distance  from  common  volume  to  receiver 

s = mean  distance  from  transmitter  to  common  volume 

g^  (r* ) = • combined  transmit/receive  antenna  gain  characteristic 
for  /.th  diversity  channel,  g^  < 1, 


§(•)  = refractive  index  spectrum 

m = unit  vector  in  direction  of  transmit  energy 
n — unit  vector  in  direction  of  scattered  energy 

U(r',t)  = 6(!r'-^|  + |r-r’l  - ct  ). 
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Some  of  the  troposcatter  communication  parameters  of  interest 
to  be  computed  from  the  cross-channel  multipath  profile  are 
received  signal  power,  2a  multipath  spread,  and  antenna  coupling 
loss.  The  received  power  for  the  ^th  diversity  channel  is  given 
by 

= I 

2ttVp 

= “TT~^  J ^ [l<(ni-n)]  d^r'  (2.23) 

^ r s 

V 

For  convenience  let  ^ = 1 represent  the  main  beam  diversity 

-1 

branch  in  an  angle  diversity  system.  Then  although  is  a 

measure  of  the  absolute  path  loss,  it  will  not  be  used  in  our 
calculations  because  it  lacks  the  empirical  factors  present  in 
NBS  and  other  path  loss  prediction  techniques.  NBS  Technical 
Note  101  [2,3]  will  be  used  for  calculation  of  P^  However, 
the  relative  path  loss  determined  from  (2,23)  will  be  a 

critical  parameter  in  the  assessment  of  the  utility  of  vertical 
or  horizontally  displaced  angle  diversity  beams.  The  signal-to- 
noise  ratio  for  a digital  system  is  defined  as  the  mean  received 
energy  per  bit  divided  the  noise  power  in  a 1 Hz  bandwidth.  The 
relationship  between  SNR  and  the  median  P in  milliwatts  is 

SNR^  (dB)  = 10  log  P^  + 1.6  - (n^  + ^B  " <2.24) 
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where 

N = noise  figure  (dB) 

F 

P = receive  antenna  feeder  loss  (dB) 

L 

R = 10  log  (r  ) , r = bit  rate  in  bits/second 
B B B 

1.6=  dB  factor  converting  median  to  mean  for  a Rayleigh 
channel  . 

For  an  analog  system  (2.24)  defines  the  SNR  in  an  RP  channel 
bandwidth  of  RgHz. 

The  2c  multipath  spread  which  determines  the  intermodula- 
tion distortion  in  an  analog  system  and  the  intersymbol  inter- 
ference and  implicit  diversity  in  a digital  system  is  defined 
by 

' / H - -h  [ 

a/  -t  o o 

Antenna  coupling  loss  results  from  an  inability  to  collect 
all  the  scatter  energy  when  high  gain  narrow  beam  antennas  are 
used.  Since  in  our  definition  of  the  combined  antenna  pattern 
g (£' ) the  maximum  gain  is  fixed  at  unity,  the  antenna  coupling 
loss  for  the  /^th  diversity  branch  is  simply 


1) 


Cl^  = J lg^(r')l^  f [lc(m-n)  ]d^r'/J  §[k(m-n)]d^r'  (2.26) 

V V 

The  refractive  index  spectrum  to  be  used  in  these  calcula- 
tions is  generally  accepted  [2.2,  2.4,  2.5]  to  have  the  form 
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§ rk  (m-n)  ] = C (z)  ] m-n|  ^ 


(2.27) 


where  C(z)  is  a height  variation  which  reflects  the  smaller  mean 
square  value  of  refractive  index  fluctuations  in  the  thinner  at- 
mosphere at  larger  height  values  z.  The  unit  vectors  m and  n 
have  the  dot  product 

m . n = cos  0 

where  0 is  the  scatter  angle  at  the  point  . The  magnitude 
of  the  vector  m - n is  then 

I m-n]  = aJ1-2  cos  6 + 1 = 2 sin  9/2  = 0. 

Thus  § takes  the  form 


$ = C(Z)  (ke)' 


(2.28) 


The  scatter  parameter  n is  normally  positive  and  thus  its 
effect  is  to  reduce  the  energy  scatter  from  points  with  larger 
values  of  0.  Typical  values  of  n range  from  unity  to  five  with 
n = 11/3  the  most  commonly  used. 


2.3  Asymptotic  Estimates 


For  the  uncorrelated  scattering  channel  with  a refractive 
index  spectrum  of  the  form  of  (2.28) ^ asymptotic  estimates  of 
multipath  spread  and  coupling  loss  can  be  obtained  in  explicit 
form  for  the  limiting  case  of  small  beamwidth.  The  main  value 
of  these  estimates  is  to  functionally  determine  the  2a  multi- 
path  spread  and  coupling  loss  dependence  on  the  effective 
earth  radius  factor  K and  the  refractive  index  spectrum  slope  n. 
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The  variability  in  2o  and  C_  can  then  be  assessed  from  knowledge 
of  the  long  term  variations  in  K and  n.  The  estimates  are  also 
useful  for  approximations  to  the  actual  value  but  for  purposes 
of  predicting  median  communication  performance,  the  detailed 
numerical  integration  of  the  scatter  volume  as  described  in  the 
next  section  is  recommended  for  determining  absolute  values. 


When  the  antenna  bearawidths,  3^=P^=^,  are  much  smaller  than 
the  minimum  scattering  angle  0^,  the  refractive  index  spectrum 
does  not  change  significantly  in  the  volume  integration.  The 
multipath  profile  and  received  power  expressions  become 


-n 


2 .3 


(t)  = j U(r',t)lg^(r')rd  r’,  B«0^  (2.29) 

V 


r s 


o 2,  2-n 
2tt  k 

p s 

I 2 2 
r s 


I Is, 


(2.30) 


V 


In  order  to  estimate  the  multipath  spread  and  coupling  loss, 
we  must  approximate  the  integrals  in  (2.29)  and  (2.30).  In 
Fig.  2.3,  four  scattering  points  are  indicated:  the  origin  "0" 

corresponding  to  the  minimum  scattering  angle  0^,  the  point  "A" 
which  is  vertically  displaced  a distance  h from  "0",  the  point 
"B"  which  is  horizontally  displaced  a distance  v;  from  "0",  and 
the  point  "C"  which  is  longitudinally  displaced  a distance  I 
from  "0",  When  the  beam  angle  y is  a quarter  beamwidth  g/4,  the 
points  A and  B correspond  to  half-power  scatter  returns  relative 
to  the  return  from  the  origin.  The  half-power  distances  are 


^1  " 8 
- <30 
'^l  ■ 8 


?«0. 


3«0 
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When  Y = P/2/  the  point  C represents  a half-power  return, 
and  the  half  power  longitudinal  distance  is 


o 


A reasonable  assumption  for  the  joint  antenna  pattern  char- 
acteristic is  a Gaussian  fall-off  in  each  direction  from  the 
origin  except  for  points  below  the  origin  from  which  no  return 
is  made.  Under  this  assumption  the  common  volume  integration  in 
(2.30)  is 


.2  „ 2 « 


2 , 2 ® 


2,2 


^ ^ j’  ^-O.eSzVhi  r ^-0.69yVwi  (*  g-0.69x Vii 


3 3 

V = 4.85  . 0.0189  Y" 

o 


(2.31) 


We  now  replace  in  (2.30)  the  free  space  loss  per  unit  distance 

4 

Pq  P^(ttD/\)  4tt,  where  D is  the  antenna  reflector  diameter  assumed 
the  same  at  each  terminal  and  then 


p „p  . 

^4 


,3  3 -n-.l  -n 

<3  0 0 f , 9 « 9 • 

o o 


Where  the  first  factor  is  the  transmit  power,  the  second  factor 
accounts  for  the  isotropic  free  space  path  loss,  the  third  facbor 
represents  the  antenna  gain  dependence,  and  the  last  factor  rep- 
resents the  volume  integration  of  the  refractive  index  inhomogeneitie.s 
Since  the  antenna  beamwidth  B is  approximately  \/D,  after  combining 
factors  one  has 


\ a p.,  f'-"  D / a,  9 « 


(2.32) 
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For  turbulent  scatter#  the  spectrum  slope  n is  11/3  which 

-2/3  -14/3 

gives  a f power  fall-off  with  frequency  and  a 6^  power 

fall-off  with  scatter  angle. 


The  received  power  dependence  indicated  in  (2.32)  includes 
the  coupling  loss  effect.  We  will  examine  this  quantity  separately 
but  first  we  estimate  the  2a  multipath  spread. 

The  surfaces  of  constant  delay  in  (2.29)  are  prolate  spheroids, 
We  can  avoid  the  difficult  integrations  by  some  simple  approxima- 
tions. The  differential  delay  for  the  points  A,  B,  and  C in 
Fig.  2.3  are  readily  computed  to  be 


At,  = d39  /4c 

8«e 

A o 

o 

ATg  = dsVsc 

3«eo 

= dge^/Sc 

9«G^ 

The  delay  shell  through  the  point  A will  contribute  close 
to  the  maximum  power  to  the  multipath  profile.  Thus  to  obtain 
an  estimate  of  the  2a  spread,  v;e  assume  a triangular  shape  for 
Q (t)  with  a base  width  of  d^6^/2c  seconds.  The  estimated  spread 
is  then 


2a  = 


2^  c 


g«9. 


(2.33) 


The  main  effect  of  the  profile  shape  assumption  is  to  fix 
the  constant  1/2./6.  For  narrow  beam  antennas  the  estimate  (2.33) 
compares  favorably  with  the  predicted  value  from  numerical  in- 
tegration procedures  described  in  Section  3.  For  the  RADC 
C-band  test  link  with  0°  elevation  angles,  we  find 

9 = 0.0094 


0 = 0.0318 

o 

2a  (estimated) 
2a  (predicted) 


55  ns. 
44  ns. 


(See  Section  3) 
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Note  that  for  small  beamwidths  the  multipath  spread  is 

★ 

independent  of  spectrum  slope  n and  inversely  proportional  to 

the  effective  earth  radius  factor  K as  6 = d/KR  . 

o o 

The  multipath  estimates  are  generally  a little  larger  than 
the  values  obtained  by  numerical  integration  of  the  common  volume 
because  in  the  integration  the  6 ^{r')  fall-off  with  height 
reduces  the  effective  multipath  spread. 


When  the  antenna  beamwidths  are  large  compared  to  the  minimum 
scattering  angle,  the  volume  integration  is  primarily  dependent 
on  the  refractive  index  spectrum,  §.  We  have 


(t) 


2tt  h P, 


2 2 
r s 


fc(2) 

V(t) 


6"^(r')d^r', 


(2.34) 


2 2 

2tt  k P 


2 2 
r s 


J C(z)  e”^(r' )d^r' , 
V 


(2.35) 


The  height  dependence  typically  falls  off  inversely  propor- 
tional to  height  and  is  not  as  a significant  factor  in  the  volume 
integration  as  the  angle  dependence.  Neglecting  this  height 
dependence  by  setting  C(z)  = C^,  a constant,  we  can  compute  the 
effective  scattering  vclume  size.  Let  0^  be  the  value  of  6 
such  that  the  scattered  power  is  reduced  by  3dB  from  the  maximum 
scattered  power,  i.e.  solves 


or 


O 


V 


o 


only  when  the  take-off  angles 


=0,  otherwise  8 =d/KR  -i-s^  + /i  . 
r o o t r 


Now  let  and  w^  be  the  height  and  width  shift  required 
from  the  point  of  minimum  scattering  angle  to  produce  a scattering 
angle  of  0^  in  the  respective  vertical  or  azimuthal  planes.  The 
geom.etry  for  the  determination  of  h^  and  is  shown  in  Fig.  2.4. 
These  values  are  easily  shown  to  be 


h 


1 


de-  de  de  , i . 


(2.36) 


(2.37) 


The  length  of  the  scattering  volume  is  on  the  order  of  the 
path  length  as  the  geometry  shov/s  that  the  longitutinal  distance 
is 


2h, 


f ^ 

"1  9,-0. 


d/2 


To  estimate  the  size  of  the  scattering  volume  assume  that  the 
vertical  and  azimuthal  fall-off  as  a function  of  departure  from 
the  point  of  minimum  scattering  is  Gaussian  distributed.  If  the 
length  of  the  scattering  volume  is  taken  as  d,  the  size  of  the 
volume  in  (2.35)  becomes 


„ , r -0.69z^/h.^,_  P -0.69yVwT^- 

Vg=dje  'IdZje  ■''Idy 

° o -oo 

= d(|  y2nw^Vl.38 

= 2.27  dh^w^  = 0.14  d^0y(2^^"-l)  (2.38) 


i 
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SIDE  VIEW 


[ 


Fig.  2.4  Scatter  Volume  Half-Power  Dimensions  | 
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(2/39) 


The  received  power  is  approximately 


P 


2 2 

2tt  k P C 
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The  asymptotic  coupling  loss  for  narrow  beam  antennas  can 
now  be  estimated  as  the  ratio  of  scattering  volumes  for  the  two 
examples  considered.  We  have 


I 


-n  3 , 
e dr* 


" . -n*' 


• \ I 


9q”  j Ig  (£' ) I r 


(2.40) 


which  shows  that  the  asymptotic  coupling  loss  varies  inversely 
as  the  cube  of  the  beamwidth  when  the  transmitter  and  receiver 
beaiTiwidths  are  constrained  to  be  the  same.  For  different  values 
of  transmitter  and  receiver  beamwidths,  the  coupling  loss  is 


C 


L 


7.5  (2^/^-l) 


p p : 

max  min 


P . <P  «6  . 

min  max  o 


If  one  of  the  two  antenna  beams  is  on  the  order  of  the 
scattering  angle,  i.e.  a wide  beamv/idth,  the  coupling  loss  at 
the  single  narrow  beam  antenna  is 


C 


L 


single 

Antenna 


P«9 

o 
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which  agrees  well  with  Gjessing's  result  T2.  5 ] derived  from 
correlation  distance  calculations  at  the  receiver  antenna.  He 
found 


single 

Antenna 


where 

A = antenna  aperature  area 

R = earth  radius 
N o 

Since  fl  = 7TrX/6A/4A/Tr  and  0 = d/R  , his  result  reduces  to 

o o 


single 

Antenna 

The  difference  in  constant  multipliers  of  1.4dB  less  loss 
for  Gjessing's  result  is  due  to  choice  of  approximation  functions. 

As  a numerical  check  on  the  coupling  loss  for  a syrometical 
link,  the  C-band  RADC  test  link  considered  in  Section  3 has 
parameters 

6 = 0.0318 

o 

B = 0.0094 

which  for  a spectrum  slope  of  n = 11/3  yields  a coupling  loss 
estimate  (2.40)  of  15.9  dB.  As  will  be  indicated  in  the  next 
section,  the  estimate  is  within  the  range  of  the  numerical  in- 
tegration prediction  technique  and  other  semi-empirical  methods 
examined. 
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SECTION  3 

TROPOSCATTER  PARAMETER  PREDICTION 


The  theoretical  basis  for  the  uncorrelated  scattering  model 
of  the  troposcatter  communication  channel  has  been  presented  in 
the  previous  section.  In  this  section  we  convert  these  theoret- 
ical concepts  into  a functional  prediction  technique  for  deter- 
mining multipath  spread#  multipath  profiles  and  the  inherent 
correlation  between  delay  cells,  and  aperture  to  medium  coupling 
loss.  The  prediction  technique  consists  of  a numerical  three 
dimensional  integration  of  the  scatter  volume.  An 'absolute 

reference  for  the  main  beam  path  loss  will  be  determined  from 

* 

computerized  NBS  calculations  and  relative  power  calculations 
for  angle  diversity  beams  will  be  derived  from  the  volume  inte- 
gration. The  volume  integration  is  not  used  for  the  main  beam 
path  loss  because  it  lacks  the  empirical  data  inputs  in  the  NBS 
calculations. 

A validation  of  the  prediction  model  will  be  accomplished 
by  comparing  multipath  predictions  with  measured  data  collected 
on  the  MDTS  field  tests  at  Rome,  N.Y.  in  the  winter  of  1976. 

After  an  analysis  is  completed  which  shows  good  agreement  between 
multipath  predictions  and  measurements,  a similar  analysis  is 
planned  to  resolve  discrepancies  between  predicted  and  measured 
aperture  - to-medium  coupling  loss  values.  Another  important 
prediction  technique  output  will  be  the  cross-channel  multipath 
profile  which  provides  the  multipath  correlation  data  required 
to  predict  average  bit  error  probability  for  adaptive  digital 
modems . 

* 

A.  G.  Longley,  et.al.;  "Prediction  of  Tropospheric  Radio  Transmission 
Loss  over  Irregular  Terrain,  A Computer  Method,  1968",  ESSA(NBS), 
Boulder,  Colorado,  July  1968. 
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3,1  Prediction  Model 


Experimental  data  on  the  delay  spread  on  a troposcatter  link 
are  typically  larger  than  predicted  by  theoretical  calculations. 
Theoretical  calculations  of  coupling  loss  also  exhibit  a consid- 
erable range  of  predicted  values.  These  calculations  are  normally 
based  on  a constant  gradient  atmosphere  as  well  as  a series  of 
approximations  to  the  propagation  path  and  the  antenna  pattern. 

The  calculations  presented  here  are  exact  for  a specific  refrac- 
tion profile  which  has  a nearly  constant  gradient  near  the  surface 
of  the  (spherical)  earth.  This  analysis  is  intended  as  the 
first  step  in  resolving  the  discrepancy  between  experimental 
and  theoretical  results  obtained  in  the  past.  The  effect  of 
terminal  equipment  filters  and  measurement  errors  will  also  be 
evaluated  in  the  evaluation  of  multipath  spread. 

3.1.1  A Model  for  the  Refraction  Profile 

In  regions  of  tropospheric  scatter  it  is  generally  acknowledged 
that  the  coefficient  of  refraction  decreases  linearly  with  height. 

We  shall  use  a power-law  model  featuring  an  approximately  linear 
decrease. 

Let  n(R)  be  the  refractive  index  as  a function  of  the  distance 
from  the  center  of  the  earth.  It  is  assumed  that  the  atmosphere  is 
spherically  layered.  The  refractive  index  will  be  assumed  to  be 
of  the  form 

n(R)  = n^*  (R^R)^. 

In  this  equation  R^  is  the  radius  of  the  earth,  n^  = n(R^)  in  the 
refractive  index  at  the  surface,  and  y is  a parameter  determined 
by  the  refractive  index  gradient  at  the  surface.  Y can  be  related 
to  the  gradient  of  the  coefficient  of  refraction  (usually  expressed 
in  N-units  per  km)  as  follows 
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AN  [N-units/km]  = - Y • 10  n /R  . 

o o 


R^  is  given  in  meters,  and  n^  can  be  taken  to  be  equal  to  1.  For 


the  standard  atmosphere  the  gradient  is  40  N-units/km  corresponding 
to  Y = 0.25.  The  propagation  path  is  naturally  curved,  but  by  a 
suitable  transformation  of  the  coordinate  system  straight-line 
propagation  can  be  used.  Define,  in  a given  great-circle  plane, 
♦-coordinates ; 


* * 1 Y 1— V 

R = R (R)  = 7^  R ^ R ^ 
1-Y  o 


d = d (distance  along  the  surface) . 


With  this  definition  the  angular  distance  between  two  points  on  the 
surface  is  transformed  from  a to 

•k 

a = a(l-Y) . 

The  effective  earth  radius  is  then  defined  by 
"eff  ' 

This  definition  is  justified  by  the  following  theorem 

Theorem;  The  propagation  path  in  the  *-coordinates  is  a straight 
line. 

Proof; 

We  refer  to  Fig 3.1  for  the  definition  of  parameters  for  the 
path.  In  the  real  coordinate  system  the  path  is  governed  by  the 
following  equations; 

(1)  R • n(R)  sin  0 (R)  = R n sin  6 (Spherical  Snell's  law). 

o o o 

6 (R)  is  the  zenith  angle  at  a variable  point  of  the  path,  and  0^ 
is  the  zenith  angle  with  which  the  path  leaves  the  surface  of  the 
earth. 


(2)  R da=  tan  8 (R)  dR 


da  is  an  infinitesimal  change  in  angular  distance  along  the  path. 


These  two  equations  determine  the  path  completely.  In  terms 
of  the  *-coordinates  these  equations  become 

(1*)  R*  sin  e(R  ) = R^^^  sin  6 ) 

ic  ic  ic  ic 

(2  ) R da  = tan  e (R  ) dR  . 

This  is  just  the  equations  governing  electromagnetic  wave  propaga- 
tion in  a homogeneous  medium.  Hence  the  rays  are  transformed  in 
straight  lines  in  this  coordinate  system.  It  is  verified  below  that 
distances  on  the  surface  are  unchanged- by  the  transformation, so  it 
follows  that  the  velocity  of  propagation  in  the  *-coordinates  is 
equal  to  the  surface  velocity  in  the  original  coordinate  system. 

This  completes  the  proof  of  the  theorem. 


The  above  resuj  •.  is  applied  to  the  troposcatter  problem  by 
considering  separately  the  path  from  a terminal  (transmitter  or 
receiver)  to  a given  point  in  the  common  volume.  This  path  is 
transformed  as  above  in  the  great-circle  plane  containing  the  ter- 
minal and  the  scattering  point.  The  height  of  the  scattering  point 
in  the  *-coordinates  is 


* 

h 


eff  L 


1-Y 


- R _ r (l-Y)  h/R  - Y(l-Y)  (h/R  ) V2I 
eff  L o o _j 

= h(l-Yh/2R  ) . 

o 
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It  is  seen  that  the  height  is  changed  only  very  little.  As  an  ex- 
ample take  h = 50,000  ft.,  and  y = 0.25  (standard  atmosphere)  then 

h*  = 49,985  ft. 

* 

The  small  difference  between  h and  h can  usually  be  ignored,  but 
correction  term  will  be  used  in  the  results  to  be  presented  later. 
The  case  of  small  h verifies  that  infinitesimal  distances  at  the 
surface  are  unchanged  by  the  transformation. 

In  the  *-coordinates  it  is  now  easy  to  find  the  transmission 
delay,  the  take-off  angle,  and  the  angle  with  horizontal  at  the 
scattering  point.  These  parameters  are  the  same  in  the  original 
coordinates.  The  angle  at  the  scattering  point  is  particularly 
important  due  to  the  drastic  scattering  angle  dependence  of  the 
turbulent  scatter  cross  section. 


In  the  following  we  derive  formulas  for  the  dimensions  of  the 
scattering  volume  and  the  delay  spread.  The  purpose  of  these  sec- 
tions is  to  document  some  of  the  formulas  used  in  the  prediction 
model.  Sections  3.1.2  and  3.1.3  can  therefore  be  skipped  at  a 
first  reading. 

3.1.2  Height  and  Size  of  the  Common  Volume 

The  extent  of  the  common  volume  in  the  azimuth  direction  will 
be  considered  last.  For  the  time  being  we  can  constrain  the  analysis 
to  the  great-circle  through  the  two  terminals.  The  coordinates  are 
transformed  in  this  plane  such  that  the  rays  are  straight  lines. 

The  transmitter  and  receiver  beams  are  characterized  by  the  take- 
off angles  of  the  upper  and  lower  rays  in  the  beam.  The  notation 
used  here  and  in  Fig. 3.2  is 


«P. 


CO 


tu 

tt 


(D 


ru 


CO 


tl 


- take-off  angle  at  transmitter  of  upper  ray. 
= take-off  angle  at  transmitter  of  lower  ray. 
= take-off  angle  at  receiver  of  upper  ray. 

= take-off  angle  at  receiver  of  lower  ray. 


The  beamwidths  are 


t tu  tl 


P = CD  •“  Cp  , 

r ru 


the  take-off  angles  of  the  center  of  the  beam  are 


''t  = <'^tu  ^ 


Y = (co  + cp  , ) /2 . 
r ru  rl 


The  altitude  h of  the  lowest  point  of  the  common  volume  is  found 
to  be  determined  by 


2h'-  h‘-  h;  + h'  - h'h*  = 
r t r t 


CP  .+  cp  CD  - CO 

'"'"i 


sin(a  + + (5^^) 


+ (1+h^)  (1+h^)  {^sin sin0  + co^^+ cp^^)  - sin 


cos^  |(a  + + C0_.^) 


where  h'  = h /R  h'  = h /R  __  and  h'  = h/R  . 
r r eff  t t eff  eff 


This  can  be  approximated  by 


h + h 

h - 

2 2R  .. 

eff 


a + CP,  , + CO  , 
tl  rl 


eff  cx  /a  , ^ 

2 2 (2  ®rt  ^tl 


^ p " p ^ 2 


1 


In  particular,  if  the  heights  of  the  receiver  and  transmitter  are 


zero, 


h ~ 


\2{<i  + <o^^+  co^^V 


a /d  . 
2 \2 


^ + CD  . 
2 ‘ r-t 


* “m) 


In  the  simple  case  where  co  , = cp  , = m,  and  h = h ~ 0,  we 

tt  r'6  t r 


have 


h - 


2 

~ ef f sina/2  • sin(Ci/2  + 2cp)  ~ ^eff 


cos 


(a  + 2rn)/2^ 


8 


- 1)- 


The  altitude  of  other  points  in  the  common  volume  are  easily  found 
by  a simple  substitution.  In  particular,  the  vertical  thickness 
of  this  volume  is  approximately,  assxaming  equal  height  terminals. 

For  a symmetric  link,  where  Y,  “ Y (equal  take-off  angle)  and 

t r 

P^=  = P we  get  simply 


A similar  calculation  of  the  length  of  the  common  volume  gives 

is  found  to  be  approximately 


The  width  of  the  volume 


L ~ d 

y 


P +P^. 

r t 
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For  a symmetric  link, 

L = L = p d/2 

y z 

~ Yj.  '*■  “)  • 

3.  1.3  Delay  Spread 

The  precise  delay  spread  will  be  determined  by  a numerical  in- 
tegration over  the  common  volume.  However,  an  evaluation  of  the 
maximum  delay  spread  within  the  volume  can  be  made  easily, given  the 
nominal  antenna  beam  widths, S. 

The  relative  delay,  as  a function  of  height  z above  the  lowest 
point  of  the  common  volum.e  is  found  to  be 


T 


Jl.  f— + 

C L2R 
o eff 


2„ 
-Z_  i 
d 


where  is  the  surface  velocity  of  propagation, 
if  z = and  the  link  is  symmetric,  we  get 


T 


A.  r Bd' 


-2R 


eff 


In  particular, 


This  represents  the  maximum  delay  spread  for  a symmetric  link.  Fig. 3. 3 
shows  the  relative  delay  t as  a function  z,and  Fig. 3. 4 shows  the 
maximum  delay  spread  as  a function  of  the  path  length  d and  the 
beamwidth  P of  a symmetric  link. 


3.1.4  Scattering  Cross  Section 


The  power  scattered  from  a small  cell  of  volume  dV  depends  on 
the  scattering  angle  and  the  spectrum  slope  parameter  n.  The  fol- 
lowing formula  for  the  power  received  from  the  volume  cell  dV  is  used: 


“ P,.  G.  (dV)  G (dV) 
r t t r 


dV  11 

n * 2*  2 

6 (dV)  R^.  Ry 


1 

h(dV) 
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The  assumptions  for  the  uncorrelated  scattering  model  and  the  deriva- 
tion of  the  received  power  equation  were  presented  in  Section  2.1. 

The  parameters  and  G are  the  gains  to  dV  from  the  transmitting  and 
t r 

receiving  antennas,  respectively  , is  the  scattering  angle, 
and  R , R^  are  the  distances  to  the  receiver  and  transmitter, 
respectively.  The  inverse  dependence  of  height  h accounts  for 
the  thinning  out  of  the  atmosphere  with  height.  This  dependence 
is  customarily  used  in  predicting  troposcatter  performances. 

3.1.5  The  Antenna  Pattern 


It  is  well-known  that  the  gain  pattern  of  a circular,  uni- 
formly illuminated  aperture  is 

G(cp)  = 2J^(C)/C 

where 


C = 2rRsin  cp/X 


I 

I 


and 


2R  = Diameter  of  the  dish 
X = wavelength, 

0 = Bessel  function  of  order  1. 
The  3 dB  beam  width  of  this  pattern  is 


In  practice,  a parabolic  dish  will  not  be  uniformly  illuminated 
due  to  the  gain  pattern  of  the  feedhorn.  Practical  antenna's 
have  a 3 dB  beamwidth  given  by 


1 
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B ~ 70° 
3dB  — D 


Rather  than  to  evaluate  the  gain  pattern  of  a dish  with  a 
specific  nonuniform  illumination,  we  have  chosen  to  use  the 
modified  pattern 

G(cp)  = 2J^(C)/C 

where 

C,  ~ 1.67  nR  sin  cp/\. 

If  the  actual  antenna  gain  need  to  be  calculated,  we  can  use  the 
expression  for  the  power  gain 


GAIN 


= 

\i.2x; 


D » X. 


The  factor  1.2  is  included  to  account  for  imperfect  illumination. 
3.1.6  The  Computer  Programs 

The  computer  programs  that  have  been  developed  under  this 
contract  to  find  the  delay  spread  on  a troposcatter  link  uses  a 
three  dimensional  integration  over  the  common  volume.  At  each 
scattering  point  the  off  boresight  angles  of  each  antenna  are 
calculated  using  exact  formulas,  and  the  gain  evaluated  as 
described  in  Section  3.1.5.  The  scattering  angle  in  the  common 
volume  is  also  found  exactly,  and  used  to  determine  the  scattered 
power  (Section  3.1.5).  The  model  also  yields  the  aperture— to— 
medium  coupling  loss  by  calculating  the  received  power  for  the 
desired  antennas  with  the  power  received  with  practically  omni- 
directional antennas.  The  result  of  these  calculations  will  be 
described  in  Sectxon  3.3.  The  application  to  the  angle— diversitv 
mode,  and  the  associated  computer  program  is  described  below. 
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3.1.7  Prediction  of  Angle  Diversity  Path  Parameters 


In  the  preceding  sections  a model  for  troposcatter  prediction 
was  described.  This  model  can  also  be  applied  to  the  angle 
diversity  problem.  The  computer  program  used  to  evaluate  the 
model  in  Section  3 has  been  modified  to  account  for  two  receiving 
beams . This  program  calculates  the  power  received  by  each  beam 
in  a given  delay  cell,  as  well  as  the  correlation  between  the 
signals  received  in  the  same  delay  cell.  This  information  makes 
it  possible  to  evaluate  the  performance  for  a receiver  using  a 
tapped  delay  line,  such  as  MDTS.  The  angle  diversity  program 
is  still  in  the  development  stage,  and  extensive  computational 
results  will  be  presented  in  a later  report.  Here  we  describe 
some  of  the  major  concerns  that  have  to  be  addressed,  as  well 
as  some  preliminary  numerical  results. 

The  two  receiving  antennas  can  be  considered  to  be  colocated, 

so  the  same  power  is  scattered  towards  each  of  them.  Let  g^^Co) 

and  g2(6)  he  the  voltage  gain  patterns  of  the  two  receiver  beams. 

Similarly,  the  off  bores ight  angles  associated  with  a small  cell 

dV  in  the  common  volume  are  denoted  9^  (dV)  and  0^  (dV) . If  the 

power  scattered  to  the  receiver  from  the  volume  cell  dV  is  P (dV) , 

s o 

the  total  power  received  on  the  first  (primary)  beam  g^(*)  is 

'’i  “ J I®!  <ei(av))l^dv 

common 

volume 

On  the  second  (angle  diversity)  beam  92^*^  received  poi\7er  is 

common 

volume 
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The  correlation  between  the  signal  received  on  the  two  beams  is 

^19  = fonn,n,on  <3ABAdv))  g*(9_(dv))  dv. 

12  V common  sc  i i 2 2 

volxime 

This  is  in  general  a complex  number,  and  it  is  important  to  realize 
that  knowledge  of  the  phases  of  the  gain  patterns  is  required  to 
evaluate  this  integral  accurately- 

In  the  results  we  have  obtained  so  far  real  gain  functions  of 
the  form  described  in  Section  3.1.5  have  been  assumed.  An  assessment 
of  the  error  increased  by  making  this  approximation  will  be  made 
at  a later  time. 

It  should  also  be  pointed  out  that  the  sidelobes  can  make  a 
considerable  contribution  to  the  correlation.  This  is  in  contrast 
to  the  single  beam  calculations  in  Section  3,1.6,  where  the  sidelobes 
were  justifiably  ignored.  This  fact  is  illustrated  in  Pig.  3.5, 
where  the  upper  and  lower  beams  are  overlapping  at  the  3 dB  point. 

The  amplitude  of  the  "cross  beam"  voltage  gain  ( (9)g*  ('^)  | ^) 
is  shown  together  with  the  two  primary  beams.  It  should  be 
clear  from  this  figure  that  the  sidelobes  of  this  "product  beam"  are 
much  more  prominent  than  those  of  the  two  primary  beams.  Hence  it  is 
more  important  to  include  the  sidelobes  in  the  cross  correlation 
calculation  than  in  the  determination  of  the  power  received  on  each 
beam. 


Due  to  the  high  complexity  of  the  computer  program,  only  a few 
cases  have  been  evaluated  at  this  time.  In  particular,  for  the  4.5 
GHz  path  in  Table  3.2  and  the  two  beams  overlapping  at  the  3 dB  point 
we  have  the  results  in  Table  3.1  below.  From  these  two  cases  it  is 
seen  that  upper  beam  loses  power  relative  to  the  lower  beam  at  the 
higher  elevation  angle.  However,  it  must  be  recognized  that  the 
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Path  Data:  4.5  GHz  Path  in  Table  3.2 


Elevation  of  Transmitting  and  | 

Lower  Receiving  Beam  i 

0° 

Effective  Earth  Radius  factor 
K 

1.3 

Spectrum  Slope  Parameter  n 


Power  Received  on  Upper  Beam, 
Relative  to  lower  beam 


-0.5  dB 


-2.6  dB 


Correlation  between  Signals  Received 
on  the  Lower  and  Upper  Beams 

0.58 

0.33 

2cr -delay  Spread,  Lower 
Beam 

46 

nsec 

93 

nsec 

2a “delay  spread.  Upper 
Beam 

67 

nsec 

118 

nsec 

power  in  the  lower  beam  is  higher  by  2-3  dB  at  the  higher  eleva- 
tion angle  (see  Fig.  3.6) . Hence  it  appears  that  elevating  the 
beam  by  a quarter  to  a half  beairPA^idth  both  increases  the  power 
and  decreases  the  correlation.  A more  extensive  example  study  of 
the  angle  diversity  problem  will  be  performed  later. 

3 . 2 Multipath  Computations 


The  prediction  model  outlined  in  Section  3.1  has  been  tested 
with  typical  data  for  which  empirical  results  have  been  obtained. 
Table  3.2  shows  the  path  configuration  and  the  two  frequencies  ana 
lyzed,  together  with  multipath  spreads  measured  during  the  MDTS 
test  program  at  the  RADC  Verona-young s town  Link  [3.1].  Several 


parameters  are  unknown  for  the  measured  paths,  such  as  the  precise 
take-off  angles  and  the  atmospheric  parameters.  The  troposcatter 
prediction  model  has  therefore  been  applied  for  several  of  those 
parameters.  The  take-off  angles  have  been  varied  to  maximize  the 
received  power,  and  the  take-off  angles  resulting  in  less  than  a 
3 dB  power  loss  have  been  found  when  possible.  In  this  way  the 
delay  spread  and  coupling  losses  can  be  found  for  both  optimally 
adjusted  antennas,  and  for  slightly  misad justed  antenna  pointing. 
The  effective  earth  radius  and  the  refractive  index  spectrum 
slope  have  been  varied  between  the  10  and  90  percentiles  found  in 
most  observations.  These  numerical  results  can  then  reliably 
predict  the  variability  in  troposcatter  propagation. 

3.2.1  Theoretical  Results 

Figure  3. 6 illustrates  the  dependence  on  take-off  angle  for 
the  900  MHz  path  in  Table  3. 2. For  full  understanding  of  these 
curves  it  should  be  added  that  the  3 dB  beamwidth  of  the  antennas 
are  approximately  2.7°.  With  a 0°  beam  elevation  half  of  the  beam 
is  blocked  by  the  earth  bulge,  so  the  received  power  is  increased 
somewhat  when  the  beams  are  slightly  elevated  , However  as  the 
beams  are  pointed  higher  the  scattering  angle  is  increased  and  this 
eventually  causes  a decrease  in  the  received  power.  The  curve  for 
the  received  power  is  quite  flat  but  the  maximum  is  seen  to  be  at 
approximately  a quarter  of  the  3 dB  beamwidth.  The  received  power 
is  reduced  3 dB  with  respect  to  its  maximum  value  at  a 1.7°  beam- 
elevation.  The  most  important  result  of  Fig.  3.6  is  the  strong 
dependence  of  delay  spread  with  the  elevation.  The  delay  spread 

varies  from  260  nsec  to  610  nsec  as  the  elevation  of  both  trans- 

. . , o o 

ma.tter  and  receiver  varies  from  0 to  1.7  . 
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Figure  3.7  shows  the  same  two  curves  for  the  4.5  GHz  path.  The 
maximum  of  the  received  power  is  even  flatter  in  this  case,  not 
falling  3 dB  below  the  peak  value  until  almost  a full  beamwidth 
elevation  (0.54°  3 dB  beamwidth) . The  delay  spread  varies  from 
44  nsec  to  110  nsec  when  the  received  power  varies  over  a 3 dB  range. 

Figure  3.8  shows  how  the  delay  spread  depends  on  the  effective 
earth  radius  factor  K and  the  spectrum  slope  of  the  refractive 
index  at  4.5  GHz.  At  this  frequency  the  beam  is  so  narrow  that 
only  a slight  dependence  on  the  spectrum  slope  is  observed.  The 
dependence  on  the  K- factor  can  be  explained  as  follows:  For  in- 

creasing K the  scattering  angle  decreases  and  hence  also  the  delay 
spread.  Clearly  the  delay  spread  also  decreases  with  increasing 
spectrum  slope,  since  this  implies  a decreasing  effective  common 
volume.  Figure  3.9  shows  the  same  results  at  900  MHz.  As  can  be 
expected,  a much  stronger  dependence  on  the  spectrum  slope  is 
found,  due  to  the  wider  beamwidth. 

Figures  3.10  and  3.11  show  the  shape  of  some  typical  multipath 
profiles.  The  characteristic  sharp  build-up  of  power  at  small  de- 
lay and  slow  fall-off  at  long  delays  is  clearly  observed  on  both 
of  these  figures. 

In  the  next  section  the  troposcatter  model  will  be  compared 
with  some  experimental  data. 

3.2.2  Comparison  with  Empirical  Results 

Troposcatter  multipath  tests  were  performed  by  Sylvania/ 

SIGNATRON  during  Jan. -March  1976.  The  results  of  these  tests  v/ill 
be  used  to  validate  the  model, but  first  we  will  describe  how  these, 
and  similar  measurements, have  been  obtained.  It  is  extremely  important 
to  understand  the  shortcomings  of  the  measurement  technique  since 
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of  Take  Off  Angles  at 


: Atmospheric  Parameters  at  900  MHz 
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DELAY  (msec) 

Fig.  3.11  Typical  Multipath  Profiles  at  4.5  GHz 


they  are  the  major  cause  for  the  discrepancies  between  theoretical 
and  empirical  results  found  in  previous  studies. 

The  probing  technique  is  illustrated  in  Fig.  3.12.  Let  the 
probing  waveform  (typically  a PN-sequence)  be  p(t)  and  let  the 
transmitted  wavefonn  be  p^(t), 

OD 

p (t)  = r h(r)  p(t-T)dT, 

^ •'c 

where  h(T)  is  the  impulse  response  of  the  transmitter  filter.  If 
the  time  varying  channel  impulse  response  is  denoted  g(t,r)  the 
received  signal  is 

CO 

r(t)  = f g(t,T)  p (t-T)dT  +n(t). 

n(t)  is  additive  white  Gaussian  noise.  The  theoretically  optimal 
receiver  is  a filter  matched  to  p^(t).  In  practice  the  received 
signal  is  correlated  with  delayed  version  of  p(t).  For  our  purpose 
we  can  consider  the  receiver  to  be  a filter  matched  to  p(t).  The 
output,  z(t),  of  this  filter  is 
z(t)  = Jr  (v  )p  (T-t+v)  dv 

= Jdv  JdT  g(v/T)  p^(v-t)p (T-t+v) 

T is  the  signal  duration.  Squaring  and  averaging  this  yields 
(ignoring  noise 

^ “ J^V2jdv2jdTQ  (v2^-V2»‘t)p^  ( Vj^-t)p  (T-t+Vj^)p^  (v2-'t)p  (T-t+V2)  • 

In  this  expression  Q(v, t)  is  the  tap  gain  correlation  function 
defined  by 

E(g  (t,Tj_)g(t+v,  T^))  = Q(v,  t^)6(t^-t2). 

It  is  assumed  that  the  channel  is  WSSUS  (Wide— Sense  Stationary 
Uncorrelated  Scatter) • In  practically  all  troposcatter  situations 
the  fading  is  so  slow  that  we  can  make  the  "frozen"  channel 
approximation 
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Q(v,t)  ~ Q (0,t)  - Q(t)  . 

This  gives 

00 

z^(t)  = J Q(t)  |p  (T-t+T)  I dT, 
o 

where  p(t)  is  the  signal  correlation 
p(t)  = Jp^(v)  p(v+T)dv. 

2 

If  1p(t)|  were  a Dirac  delta  function  we  would  get 
(t)  = Q (t  - T)  . 

Probing  signals  are  designed  to  approximate  such  a delta  function, 
so  we  tabe  as  the  estimate  of  Q(t), 

Q (t)  - z“(t+T). 


Usually  the  delay  spread  is  characterized  by  the  rms  spread  of 
Q (t)  » defined  by 


T^Q(T)dT/  jQ(T)dT  - [jTQ(r)dT/  J Q(T)dT]^. 


With  this  definition  the  m.easured  delay  spread  is 


I 1 2 

The  second  terms  represents  the  error  due  to  the  fact  that  !p| 
is  not  a delta  function.  This  also  includes  the  effect  of  the 
transmitter  filter. 


In  the  most  common  situation  a PN  sequence  is  transmitted, 
and  P (t)  can  be  assumed  to  be 


p (t)  = 


r/T, 


1t|  < T^ 
It!  > T 
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T is  the  period  of  the  pulses.  If  the  filter  can  be  ignored 
o 

and  an  essentially  continuous  sampling  is  assumed  we  have 

For  the  measurements  used  in  this  report  we  have 

T = 100  nsec, 
o 

It  follows  that  the  minimum  2c  multipath  spread  is 


^ “ vio  ’’o  ' 

In  practice  the  transmitter  filter  can  often  cause  long 
tails  on  p (t) , and  hence  substantially  increase  the  multipath 
spread  inherent  in  the  measurement  system.  As  a simple  example, 
Table  3.3  shows  the  effect  of  a single-pole  filter. 

Table  3.3 

2o  Delay  Spread  of  Ideal  Channel  v/ith  a Single  Pole  Filter 


2a|Q| 2/nsec 


The  effect  of  the  additive  noise  is  to  add  a constant  to  the 

A 

estimated  delay  power  spectrum  Q(t).  This  means  in  principle  in- 
finite delay  spread.  This  startling  fact  is  somewhat  alleviated 
by  realizing  that  only  a finite  number  of  samples  of  Q(t)  are  used 
in  practice,  but  it  clearly  underscores  the  importance  of  eliminat- 
ing the  noise.  Also  it  must  be  pointed  out  that  the  RAKE  adds  a 

noise  dither  signal  at  the  A/D  input  which  produces  a -15  dB  in- 
* 

put  SNR.  The  resulting  noise  of  the  output  of  the  integrator  is 

* The  output  SNR  is  -15  dB  + 10  log  T /T  where  T = integration 
time,T  = symbol  length. 
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the  main  contributor  to  the  above  mentioned  additive  noise.  In 
many  of  the  measurements  the  RAKE  tap  value  due  to  this  noise 
alone  contribute  significantly  to  the  measured  spread.  Since  only 
10  taps  were  used, the  error  is  largest  when  the  actual  delay  spread 
is  relatively  small  (say  200  nsec  or  less) . For  this  reason  it  is 
mainly  the  4.5  GHz  channel  that  is  affected,  while  the  measurements 
at  900  MHz  are  reasonably  accurate  (although  the  tail  end  often  has 
been  cut  off) . The  results  of  the  measurements  were  summarized  in 
Table  3.1,  where  the  effects  of  noise  and  measurement  bias  have 
been  corrected  for. 

A 

The  fact  that  the  prober  only  samples  Q (t)  at  the  rate  of 
the  transmitted  PN-sequence  means  that  the  assumption  of  essen- 
tially continuous  sampling  is  only  approximately  satisfied  at 
long  delay  spreads.  At  smaller  delay  spreads  variation  in  the 
sampling  times  can  cause  a large  variability  in  the  estimated 

A 

spread.  Suppose  Q(t)  is  sampled  at  the  instants  = iT^.  The 
estimated  rm.s  spread  is  calculated  by 


i 

‘ a 


A » A 

Q = Q (iT^) 


(r  a,  ^ a,  - (.  ia,)^)/  (.  aj^ 


In  general  this  does  not  relate  simply  to  a q » so  let  us  consider 


a few  simple  examples. 
Example  1; 


Q('J')  = 6 (t-t  ),  o < T < T 
o o o 


In  this  extreme  case  we  clearly  have 


- Q = 0- 
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J 

...  1 


However.  Q.  is 
1 

CO 

Qi  ■=  J Q (t)  lQ(iT^  + t)  1^  <3t 

0 

= lQ(iT„  + T^)|^ 

p i < 0 

1 (1-T  /T  ) ^ i = 0 

; O O 

i = l 

t o o 

’0  i > 1 

Hence  the  estimated  delay  spread  is  determined  by 

2-  _ n,2  ^0^1 

^ Q " (Q^Q.  ) 2 • 

o i 

If  Tq  = 0 or  Tq  we  get  the  correct  result,  0^=0, 

However  if  t = T /2  we  get 

o o ^ 

Hence,  if  T = 100  nsec  there  can  be  an  error  in  the  2a  - delay 
o 

spread  by  the  same  amount,  100  nsec.  This  is  larger  than  the 
63  nsec  for  continuous  sampling,  but,  in  contrast,  it  is  also 
possible  to  get  the  exact  result. 

It  can  be  concluded  from  this  example  that  the  sparse  sam- 

A 

pling  of  Q(t)  can  be  the  source  of  much  of  the  variation  in  the 
observed  delay  spreads  when  most  of  the  received  power  falls  in 
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a delay  interval  less  than  the  tap  spacing  Tq  • We  will  consider 
one  more  example  to  substantiate  this. 

Example  2: 


Q (t) 


( 0 
I 


T < 

O 


< T < 
O — 


T + T 
o o 


T + T < T 
o o — 


2 2 

The  corresponding  delay  spread  is  a = T^/12,  or  2a  = T^/a/3. 

For  simplicity  we  consider  only  the  cases  = 0 and  T^/2.  If 

T = 0 we  have 
o 


/* 


j 0 i < 0,  i > 1 

I T /3  i = 0,  1 
V o 


Hence 


■r^/4 

o 


The  2a  spread  is  therefore  too  large,  by  a factor  of  */5,  when 


If  T = T /2  we  have 
o o 


A 


0 

T /24 
o 

7T  /12 
o 


i < 0,  i > 2 

i = 0,  2 
i = 1 
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This  gives  an  estimated  variance  of 


Ga  j~T  = T /2 
a L o o _ 


1 = 


/8 


The  2a  spread  is  now  too  large  by  a factor  of  V2  . It  is  interest- 
ing to  note  that  if  the  above  estimates  are  corrected  as  for  con- 
tinuous sampling  we  get 

Ga  |t  =0;  corrected!  = 3tV20  = vx  • -T 
a L o _!  o'  12  5 


and 

T^ 

to  ‘ corrected  1 = T^/40  = -jf  • ^0  . 

-The  average  of  these  delay  variances  is  almost  equal  to  the  correct 
value. 

From  these  two  examples  we  conclude  that  the  correction  derived 

for  continuous  sampling  can  be  used  for  both  long  and  short  delay 

spreads,  but  that  it  must  be  interpreted  as  an  average  over  sampling 

times  at  short  delay  spreads,  with  considerably  variation  possible 

around  this  value.  In  general,  delay  spread  shorter  than  T (100 

o 

nsec)  cannot  be  expected  to  be  accurate. 

With  the  above  discussed  measurement  errors  in  mind  we  can 
now  compare  the  empirical  and  theoretical  results.  The  results 
of  the  measurements  were  tabulated  in  Table  3. 2.  At  900  MHz  the 
median  delay  spread  is  290  nsec,  if  the  error  due  to  noise  and 
sampling  are  ignored.  This  is  a reasonable  approximation  at  this 
frequency  but  an  uncertainty  of  up  to  + 50  nsec  should  be  allowed 
for.  This  measured  delay  spread  can  be  compared  with  the  theoretical  re- 
sult of  260  nsec  with  antenna  pointed  at  the  horizon  and  K=  1.3, 

N=3  (Fig. 3.6)  .If  the  antennas  are  pointed  to  maximize  the  power 
the  delay  spread  is  seen  to  be  approximately  340  nsec  with  the 
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same  atmospheric  parameters.  FromFig.  3.9  it  is  seen  that  290 
nsec  spread  is  achieved  with  0.7°  take  off  angles,  K = 1.3,  and 
N - 14/3. 

The  correspondence  between  theoretical  and  empirical  results 
at  900  MHz  is  well  within  uncertainties  due  to  unknov/n  atmospheric 
parameters,  pointing  angles,  noise,  and  time  synchronization. 

At  4.5  GHz  the  delay  spread  is  so  small  that  it  is  necessary 
to  subtract  an  estimate  of  the  noise  from  Q(t),  and  to  compensate 
for  errors  in  the  correlation  process.  As  justified  above,  the 
correction  corresponding  to  continuous  sampling  has  been  used. 

The  results  (corrected)  were  symmarized  in  Table  3.1,  showing  a 
median  delay  spread  of  118  nsec.  In  contrast,  the  median  without 
correction  is  approximately  145  nsec.  The  theoretical  value  when 
the  beams  are  pointed  at  the  horizon  is  only  44  nsec.  (Fig.  3.6) 
but  with  the  same  received  power  the  beams  could  be  pointed  0.47° 
above  the  horizon,  corresponding  to  a delay  spread  of  108  nsec. 

It  is  in  fact  likely  that  one  or  both  of  the  antennas  are  pointed 
slightly  upwards  to  avoid  local  obstructions.  The  measured  delay 
spread  can  easily  be  explained  this  way,  but  it  is  perhaps  more 
likely  that  the  delay  spreads  estimated  from  the  measurements 
are  still  too  large.  The  large  uncertainties  associated  with 
timing  errors,  receiver  noise  and  dither , pointing  angles,  and 
atmospheric  noise  exnlains  easily  the  small  differences  (50-100 
nsec)  between  the  empirical  and  theoretical  data. 

3.2.3  Summary,  Multipath  Computations 

The  theoretical  model  developed  in  Section  3.1  gives  a rea- 
sonable agreement  with  measured  results.  As  can  be  expected  the 
differences  a^re  most  pronounced  at  the  high  frequency  (narrow 
beam)  where  the  results  are  much  more  sensitive  to  given  parameter 
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variations.  The  agreement  is  well  within  the  uncertainties,  and 
much  better  than  reported  in  previous  studies  [3.2]. 

3 . 3 Coupling  Loss  Computations 

The  aperture-to-medium  coupling  loss  on  a troposcatter  link 
is  defined  as  the  loss  incurred  by  not  using  the  total  available 
scatter  potential.  In  other  words,  omnidirectional  antennas  have 
no  coupling  loss,  v;hile*the  gain  of  two  narrow  beam  antennas  is 
offset  by  the  coupling.  It  will  be  shown  that  from  elementary 
considerations  the  total  gain  increases  asymptotically  for  very 
narrow  beams  by  3 dB  each  time  the  beamwidths  are  halved.  This  cor- 
responds to  a 9 dB  aperture-to-medium  additional  coupling  loss  in- 
curred when  the  beamv;idths  are  halved. 

The  troposcatter  prediction  model  developed  herein  can  be 
used  to  predict  the  coupling  loss.  While  the  validity  of  the  model 
for  delay  spread  calculations  was  verified  in  Section  3.2,  we  com- 
pare here  the  predicted  values  of  the  aperture-to-medium  coupling 
loss  with  several  empirical  results. 

3.3.1  Predicted  Aperture-to-Medium  Coupling  Loss 

The  coupling  loss  is  calculated  simply  by  computing  the  power 
received  with  a wide,  essentially  omnidirectional  beam  and  the 
desired  narrow  beamwidths.  The  ratio  of  the  received  powers  (elim- 
inating the  antenna  gains)  is  the  coupling  loss.. 

In  Fig.  3.13  the  computed  coupling  loss  is  shown  for  the  4.5 
GHz  link  in  Table  3.1.  As  the  two  atmospheric  parameter  N,K, varies, 
the  coupling  loss  varies  between  10  and  16  dB.  This  can  be  com- 
pared with  some  known  results.  The  nomogram  (by  Collins  Radio) 
in  Figs.  13-22  in  Panter  [3.3]  yields  the  estimate  11.5  dB.  The 
curves  (by  Hartman  and  Wilkerson,  1959)  in  Figs.  12-15  also  in 
Panter  tells  us  that  a 10  dB  coupling  loss  can  be  expected. 


3-37 


These  numbers  are  all  in  the  low  end  of  the  range  of  the  coupling 
loss  calculated  from  our  model.  Hov/ever,  it  must  be  realized 
that  all  these  nomograms  and  curves  are  constructed  from  empirical 
data,  where  the  antennas  may  well  be  pointed  higher  than  the  0° 
elevation  angle  assumed  in  our  calculations.  This  uncertainty  in 
the  reported  measurements  places  the  actual  coupling  loss  for  the 
link  configuration  use  in  Fig. 3.13  somewhere  in  the  range  of  10-15 
dB,  in  better  agreement  with  the  predicted  results. 

The  results  in  Fig.  3. 13  are  for  the  4.5  GHz  link.  Figure  3.14 
shows  the  aperture-to-medium  coupling  loss  as  a function  of  fre- 
quency. The  calculations 'have  been  performed  for  two  relatively 
extreme  sets  of  atmospheric  parameters.  In  combination  with 
Fig.  3 . 13  it  is  possible  to  find  the  coupling  loss  by  extrapolation 
for  almost  my  situation.  Also  shown  in  Fig. 3 . 14  are  some  of  the 
known  results  referenced  in  connection  with  Fig.  3.13. 

It  should  be  emphasized  that  these  predicted  values  of  coupling 
loss  represent  median  quantities.  The  long  term  variability  of 
coupling  loss  can  be  computed  from  Figs. 3. 13  and  3.14  given  prob- 
ability distribution  information  on  the  effective  earth  radius 
factor  K and  the  refractive  index  spectrum  slope  N.  The  knowledge 
of  antenna  take-off  angles  and  atmospheric  parameter  variations  is 
essential  if  the  effects  of  aperture-to-medium  coupling  loss  are 
to  be  adequately  accounted  for. 

3.3.2  Specific  Considerations 

Coupling  loss  effects  have  special  significance  for  transmis- 
sion systems  where  radio  carrier  frequency  changes  are  anticipated 
because  the  increased  loss  at  higher  frequencies  may  preclude  opera- 
tion with  adequate  system  margin.  We  examine  here  some  of  the 
salient  features  associated  with  a frequency  change  (doubling  of 
the  frequency  is  used  as  an  example) , provide  a simplistic  but 
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asymptotically  accurate  interpretation  of  coupling  loss,  and 
finally  we  discuss  the  relative  advantages  of  multiple  feedhorns 
(angle  diversity)  vs.  multiple  antennas  (space  diversity)  as  a 
means  of  alleviating  the  coupling  loss  problem. 

If  the  operating  frequency  of  an  existing  tropo  link  is 
doubled  there  are  three  factors  affecting  received  power: 

(a)  Free-space  antenna  gain  at  each  terminal  increases  by 
6 dB  for  a total  increase  of  12  dB. 

(b)  The  tropo  path  loss  between  isotropic  radiators  in- 
creases by  9 dB.  (See,  for  example,  P.F.Panter,  Com- 
munication Systems  Design, pp. 384-397 , McGrav;-Hill,  1972. ) 

Of  this  9 dB,  6 dB  is  the  normal  free-space  isotropic 
aperture  area  loss,  and  3 dB  is  the  scatter  anomaly. 

(c)  The  aperture- to-medium  coupling-loss  increases. 

There  is  thus  a 3 dB  gain  excluding  the  change  in  coupling-loss. 

Using  the  nomographs  in  Panter  (op.cit.)  the  change  in  coupling 
loss  can  be  evaluated  for  a few  simple  cases,  listed  in  the  following 
table. 


TABLE  3.4 

Increase  in  Coupling  loss  from 
Doubling  the  operating  frequency 


14  foot  dish 
1.00  at  4.4GHz 
2.0°  at  2.2GHz 

28  foot  dish 
0.5°  at  4.4GHz 
1.0°  at  2.2GHz 

Path  length  100  nmi 
Scatter  angle  1.25° 

1 dB  at  2.2GHz  1 3 dB  at  2.2GHz 
3 dB  at  4.4GHz  ‘ 10  dB  at  4.4  GHz 
Change  = 2 dB  ' Change  = 7 dB 
Net  gain:  1 dB  ; Net  Loss:  -4  dB 

Path  length  200  nmi 
Scatter  angle  2.5° 

? 

2 dB  at  2.2GHz  ‘ 7 dB  at  2.2GHz 
7 dB  at  4.4GHz  , 16  dB  at  4.4GHz 
Change  = 5 dB  : Change  = 9 dB 
Net  Loss: -2  dB  ^ Net  Loss  = 6 dB 
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The  mechanism  of  coupling  loss  can  be  understood  somewhat 
better  by  looking  at  the  shape  of  the  useful  scattering  region.  This 
region  is  bounded  below  by  the  tangent  planes  to  the  4/3  earth 
at  transmitter  and  receiver,  which  determine  the  lower  limits  of 
common  visibility.  It  is  (roughly  speaking)  bounded  above  by  a 
curved  surface  where  the  scattering  angle  is,  say,  twice  the  minimum 
value.  This  curved  surface  is  generated  by  passing  a circle  of  a 
certain  specific  radius  through  the  transmitter  and  receiver  and 
rotating  it  while  holding  those  two  points  fixed.  The  resultant 
scattering  region  looks  something  like  that  shown  in  Figure  3.15  in 
plan  and  elevation.  It  should  be  noted  that  the  scattering  ef- 
ficiency is  greatest  on  the  shortest  path  through  the  point  marked 
A on  the  figure,  and  falls  off  several  dB  toward  the  edges. 

From  the  figure,  we  see  that  as  the  dimensions  of  the  trans- 
mitted beam  get  very  narrow,  the  actual  scattering  volume  defined 
by  the  transmitter  becomes  just  proportional  to  the  product  of 
the  azimuth  and  elevation  beamwidths  . Thus  if  we  halve  the  two 
beamwidths,  we  quadruple  the  incident  power  density  but  quarter  the 
illuminated  volume  thus  scattering'  the  same  total  power  out  of  the 
beam.  Now  if  the  receiver  has  the  same  size  antenna,  there  is 
essentially  no  loss  of  power  due  to  its  azimuth  beamwidth,  which 
is  seeing  the  same  width  of  scattering  region  that  the  transmitter 
illuminates.  However,  if  the  elevation  beamwidth  of  the  receiver 
is  halved,  only  half  as  much  length  of  the  illuminated  region  is 
seen.  The  consequence  is  that  a total  increase  of  free  space  gain 
of  12  dB  results  in  only  a 3 dB  increase  in  received  power,  so  the 
coupling  loss  asymptotically  increases  by  9 dB  for  every  6 dB  in- 
crease in  dish  gain.  Another  way  of  viewing  this  is  to  recognize 
that  the  common  scattering  volume  is  proportional  to  the  cube  of 
its  linear  dimension,  or  inversely  proportional  to  the  cvibe  of  dish 
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Fig. 3. 15  Sideview  and  Top  View  of  Scattering  Region 


diameter  in  wavelengths,  so  that  since  total  gain  is  proportional 
to  fourth  power  of  dish  diameter,  the  net  gain  is  only  increasing 
linearly  with  the  dish  diameter  in  wavelengths.  With  this  way  of 
looking  at  the  loss  phenomenon,  we  can  view  the  total  dB  coupling 
loss  as  consisting  of  three  components. 

(a)  1/3  of  the  dB  loss  from  the  vertical  beamwidth  of 
the  transmitter 

(b)  1/3  of  the  d3  loss  from  the  vertical  beamwidth  of  the 
receiver 

(c)  1/3  of  the  dB  loss  from  the  narrower  of  the  transmitter 
and  receiver  azimuth  beamwidths. 

The  preceding  argument  shows  that  roughly  1/3  of  the  dB  coupling 
loss  can  be  regained  by  using  dual  vertical  angle  diversity  at  the 
receiver.  This  would  result  in  the  loss  values  of  the  following 
table . 

Table  3.5 

Median  Loss  Incurred  from  Doubling  the  Operating 
Frequency  and  Using  Vertical  Angle  Diversity 


14  foot  dish 

28  foot  dish 

100  nmi 

Gain;  2 dB 

Loss  2 dB 

200  nmi 

Loss;  0 dB 

Loss  3 dB 

The  median  loss  is  computed  for  an  idealized  vertical  angle 
diversity  system  with  equal  but  uncorrelated  energy  in  the  elevated 
beam.  For  such  a system  there  would  also  be  a diversity  gain  advant 
age.  This  diversity  gain  advantage  is  reflected  in  Fig.  3.15  where 
ideal  CPSK  at  fourth  and  eight  order  diversity  is  compared.  At 
10  ^ bit  error  probability  there  is  a total  improvement  of  approxi- 
mately 5 dB,  3 dB  of  which  is  the  median  loss  improvement  and  2 dB 
is  the  diversity  gain  advantage.  Note  that  no  additional  power 
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must  be  transmitted  in  an  angle  diversity  system  in  order  to 

.-5 


realize  this  5 dB  improvement  at  10  BER.  The  improvement  re- 
sults from  an  interception  of  additional  transmitter  beam  energy 
not  used  by  the  quad  diversity  system.  Of  course,  realistic  sys- 
tems will  not  realize  perfect  decorrelation  and  exactly  equal 
energy  in  the  elevated  beam  to  which  also  must  be  included  the 
effects  of  multipath  (implicit  diversity)  and  implementation  losses. 
Thus  a comprehensive  performance  evaluation  must  be  undertaken  in 
order  to  assess  the  full  impact  of  frequency  changes  on  existing 
systems . 


In  any  case,  for  any  fixed  order  of  diversity  and  receiver 
complexity  at  the  nev/  frequency,  there  is  probably  no  technique 
which  is  significantly  more  efficient  than  vertical  angle  diversity! 


Consider,  for  example,  the  case  of  a 140  nmi  path  using  20  foot 

.o 


aishes  at  4.4  GHz.  These  values  yield  1.75  scattering  angle  and 
o 


0.7  antenna  beamwidth  which  in  turn  results  in  an  8 dB  coupling 
loss.  With  vertical  angle  diversity,  this  could  be  reduced  to  a 
median  effective  loss  of  5 dB,  after  diversity  combining.  We  might 
consider  as  an  alternative  approach  dual  space  diversity  using  2 
dishes  with  14  foot  diameters  in  place  of  every  20  foot  dish;  the 
total  aperture  area  is  then  the  same  as  the  angle  diversity  setup. 
The  coupling  loss  with  this  setup  is  reduced  to  4h  dB.  However, 
the  total  antenna  gain  is  decreased  by  6 dB  yielding  a total  ef- 
fective loss  per  antenna  of  10%  dB  which  only  reduces  to  7%  dB  after 
diversity  combining.  Thus  the  vertical  angle  diversity  is  2%  dB 
better  than  the  equivalent  complexity  and  aperture  space  diversity 
setup. 
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ADAPTIVE  ANTENNA  CONTROL  TECHNIQUES  i 

I 

The  Defense  Communication  System  utili?,es  transhorizon  radio 
communication  links  as  an  important  component  of  its  strategic 
communication  network.  The  propagation  m.ode  for  these  transhorizon 
link  is  either  forv;ard  scatter  or  diffraction.  The  performance  of 
both  of  these  modes  is  sxabject  to  variations  in  angle  of  arrival  ; 

due  to  changes  in  the  effective  earth  radius,  ducting  phenomena,  1 

Or  shifts  in  the  effective  common  antenna  volume.  We  summarize  I 

I 

here  three  major  approaches  for  the  compensation  of  this  effect.  j 

One  of  the  approaches,  angle  diversity,  not  only  has  the  potential 
to  compensate  for  angle  of  arrival  variation  but  can  utilize  the 
angle  of  arrival  decorrelation  to  provide  a performance  improve- 
ment due  to  increased  diversity.  This  section  concludes  with  a 
brief  review  of  previous  angle  of  arrival  experimental  programs.  | 

I 

Mechanical  Steering  1 

! 

V 

The  most  obvious  method  of  correcting  for  long  term  variations  , 

in  the  angle  of  arrival  of  the  received  signal  is  to  move  either  ! 

the  antenna  or  the  feedhorn  structure.  Control  of  the  inechanical  j 

movement  can  be  effected  by  adaptation  algorithms  which  maximize  ji 

median  received  power  or  some  monotonically  related  function  such  j 

as  signal-to-noise  ratio.  Two  types  of  control  are  possible,  < 

single  ended  (open  loop)  or  double  ended  (closed  loop) . For  single  ! 

ended  control,  each  antenna  pointing  angle  is  adjusted  independently 
to  maximize  its  median  received  power.  Adjustments  must  be  in- 
frequent, on  the  order  of  minutes  in  order  to  allow  time  for  measure- 
ment of  median  power.  The  reciprocal  nature  of  the  transmit/receive 
paths  in  a typical  transhorizon  radio  link  make  single  ended  opera- 
tion without  feedback  control  possible.  When  reciprocity  does  not 
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hold,  for  example  a single  transmit  horn  and  multiple  receive 
horms,  closed  loop  control  can  be  effected  by  using  feedback  in- 
formation transmitted  over  an  order  wire  channel.  The  control 
must  consider  the  duplex  nature  of  the  transhorizon  radio  system. 
One  approach  is  to  adjust  the  antenna  pointing  angle  one  at  a 
time  and  always  in  a direction  which  jointly  maximizes  the  re- 
ceived signal  power  at  both  ends.  In  most  applications  only  ad- 
justments in  the  vertical  plane  are  required  as  earth  bulging  or 
ducting  phenomenon  produces  greater  variation  along  the  vertical 
axis. 


Because  of  the  difficulty  of  moving  large  parabolic  antennas, 

an  attractive  alternative  for  mechanical  steering  is  the  movement 

of  the  feedhorn  which  illiminates  the  parabolic  surface.  As  can 

* 

be  seen  from  Fig.  4.1  the  loss  associated  with  an  off-axis  feed 
is  very  small  for  large  pointing  angle  variations. 

A major  advantage  of  mechanical  steering  is  the  capability 
to  point  the  antenna  at  any  feasible  signal  arrival  direction  with 
little  or  no  beam  degradation  from  bores ight  conditions.  Elec- 
tronic steering  to  be  discussed  next  is  limited  by  power  levels 
and  physical  considerations  which  under  certain  conditions  pre- 
clude antenna  pointing  at  some  angles  without  beam  degradation. 

4.2  Electronic  Steering 

Electronic  beam  steering  of  a high  gain  aperture  such  as  a 
parabolic  antenna  can  be  accomplished  by  gain  and  phase  adjustments 


Data  from  Microwave  Antennas  and  Techniques,  M.I.T.  Lab. 
Series  No.  12,  edited  by  S.  Silver,  McGraw  Hill  Book  Company, 
1949. 
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of  multiple  excitations  which  are  used  to  illuminate  the  aperture. 

As  an  example  consider  two  feedhorns  arranged  in  a vertical  stack 
as  in  Fig.  4.2.  Consider  the  top  feedhorn  to  be  directed  at  the 
parabola  center  to  produce  a beam  in  the  bores ight  direction.  The 
lower  horn  generates  an  elevated  beam  raised  a degrees  above  the 
boresight  direction.  If  H is  the  horn  size  and  F is  the  focal 
length  of  the  antenna  system,  the  angle  a is  given  by 

a=  (BEAM  FACTOR)  tan"^  (H/F).  (4.1) 

The  beam  factor  accounts  for  non  uniform  illxminatiori  of  the 
aperture  and  spill-over  beyond  the  aperture.  Parabolic  antennas 
for  transhorizon  radio  systems  have  beam  factors  on  the  order  of 
0.8.  Amplitude  and  phase  adjustments  of  the  feedhorn  signals  form 
an  antenna  beam  at  any  angle  between  0 and  c from  boresight. 

If  the  beam  separation  angle  a is  too  large  in  an  electronic 
steering  application  a significant  reduction  in  antenna  gain  can  re- 
sult at  angles  near  a/2.  The  gain  in  dB  as  a function  of  angle  p 
from  boresight  of  a parabolic  antenna  can  be  approximated 

Gp  = 20  LOG  [2J^  (3.2p/n)  / (3.2  3/n)]  = -3  (2p/n)^  dB  (4.2) 

where  Q is  the  half-power  beamwidth.  Say  we  wish  to  aim  the  anten- 
na at  the  intermediate  point,  a/2.  This  is  accomplished  by  apply- 
ing equal  amplitude  and  phase  weights  to  the  two  feedhorn  signals. 
If  single  horn  excitation  at  boresight  corresponds  to  unity  volt- 
age, then  after  splitting  the  power  equally  between  feedhorns, 
each  beam  at  a/2  contributes 0.707x l0“^  volts.  The  an- 

tenna gain  relative  to  single  horn  excitation  as  a function  of  a is 

G^  = 3 - 3(a/n)^  dB  (4.3) 
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Fig.  4.2  Multiple  Feedhorn  Excitation  of  a Parabolic  Antenna 


If  the  feedhorns  can  be  constructed  such  that  the  3 dB  beam- 
v/idths  overlap,  no  loss  in  antenna  gain  results  from  pointing  the 
antenna  midway  between  the  beams.  Physical  limitation, however , 
restrict  the  minimum  value  of  a and  may  result  in  some  less  at 
the  midv/ay  pointing  angle.  For  example  a horn  not  designed  for 
an  electronic  steering  application  has  a typical  size  equal  to  tv/ice 
the  product  of  focal  length  (F)  and  wavelength  divided  by  reflector 
diameter  (D) . Allowing  for  flange  clearance,  a practical  engineering 
choice  would  be 

H = 2.3  (F/d)X. 

Larger  horn  sizes  result  in  higher  sidelobe  levels  due  to 
phase  changes  at  the  parabola  extreme ties  and  smaller  horn  sizes 
result  in  some  reduced  antenna  efficiency  due  to  energy  spillover  at 
the  parabola  edges.  Using  (4.1)  we  can  calculate  the  angle  be- 
tween beams  for  a multiple  feedhorn  application  with  this  horn  size, 

a = 0.8  tan“^  ?.t.3(F/D)X.  ^ (2.3)X/D. 

F 

Since  the  antenna  beamwidth  In  degrees  is  approximately  60X/D, 
we  have 


^ i ^ ^ ^ ^ (4.5) 

60 

However  for  electronic  steering  this  represents  a relative 
gain  of 

G = 3 - 3 (1.76)^  = - 6.3  dB  (4.6) 

r 

which  is  excessive  for  any  practical  transhorizon  radio  system. 

We  conclude  from  this  simple  analysis  of  multiple  feedhorn  excita- 
tion that  more  sophisticated  methods  of  deriving  displaced  beams 
are  required  if  electronic  beam  steering  is  to  be  a viable  approacii. 
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Some  of  the  alternatives  for  reducing  the  beam  displacement 
for  electronic  steering  include  ^ 

• sacrifice  of  antenna  efficiency  in  exchange  for  reduced 
horn  size, 

• use  of  multimode  excitation  within  a single  horn  to 
produce  displaced  beams, 

• use  of  a dielectric  lens  in  conjunction  with  vertically 
stacked  feedhorns. 

The  reduction  in  horn  size  necessary  to  move  the  beam  close 
enough  such  that  only  a 1 dB  degradation  in  the  combined  beam 
at  the  midway  pointing  angle  can  be  calculated  from  (4.1)  and  (4.3) . 
For  1 dB  loss  wa  must  have 

- 1 = 3 - 3(a/fi)^ 
or 

a/a  = vV2  . (4.7) 

Then  from  (4.1) 

H _2i_  _ ^ _ J3(60)  \ 

F " 0.8  “ 1.6  ^ ” 1.6(57.3)  D 

H = 1.13(F/D)  X,  (4.8) 

which  is  a 49%  reduction  in  horn  size  from  the  engineering  design 
example  described. 

An  engineering  design  of  a reduced  feedhorn  aperture  is  under 
development  by  R.  F.  Systems,  Inc.,  subcontractor  to  SIGNATRON.  The 
preliminary  results  of  the  design  show  that  the  loss  in  efficiency 
due  to  beam  spillover  balances  well  with  improved  edge  intensity 
such  that  a significant  decrease  in  horn  size  is  possible  with- 
out concomitant  decrease  in  antenna  efficiency.  With  a square 
horn  size  of  0.58  \ and  a beam  factor  of  0.9,  calculations  show  only 
a 0.05  dB  loss  at  boresight.  Since  most  troposcatter  systems  have 
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an  F/d  ratio  of  approximately  0.5,  this  horn  size  will  facilitate 
electronic  steering  with  less  than  a 1 <3B  degradation.  The  beam 
separation  angle  for  a horn  size  of  0.58  X and  a F/D  of  0 . 5 would  be 


a 


. _ . -1  0.58X  . 0.8(0.58)57.3 

~T.1E  = " (0.5)  ('60) ■ 


Cl 


a = n, 

i.e.,  the  half  power  points  of  the  adjacent  beams  would  intersect. 
Because  the  feedhorn  design  developed  by  RF  Systems  successfully 
overlaps  the  beams  at  the  3 dB  power  points  with  very  small  ef- 
ficiency loss  (~0.05  dB)  , this  ctpproach  is  an  obvious  choice  for 
either  electronic  steering  or  angle  diversity.  The  other  methods 
for  small  beam  separation  are  discussed  in  subsequent  paragraphs 
for  completeness. 

Monopulse  radars  use  a vertically  polarized  and 

mode  in  a single  horn  to  provide  azimuthal  angle  of  arrival  measxxre- 
ments.  Figure  4.3  illustrates  how  horizontally  polarized  TE  modes 
could  be  used  to  produce  two  independent  displaced  vertical  beams. 

The  addition  of  a TE^^  and  TE2Q  in  phase  produces  a vertical  squint 
in  one  direction  and  the  addition  of  a and  TE^q  out  of  phase 

produces  a vertical  squint  in  the  other  direction.  A single  feedhorn 
at  the  parabola  focus  can  support  both  electromagnetic  modes.  The 
major  limitation  of  this  approach  is  the  requirement  for  trans- 
horizon systems  to  operate  with  dual  polarized  horns.  In  a qua- 
druple diversity  forward  scatter  system,  the  typical  transmit  and 
receive  configuration  is  as  shown  in  Fig.  4,4,  Note  that  the  feed- 
horn must  receive  both  vertical  and  horizontally  polarized  waves . 
The  use  of  dual  polarized  feedhorns  allows  for  more  convenient 
signal  separation  and  is  standard  practice.  Thus  although  the 
horizontally  polarized  wave  can  be  used  in  a multimode  horn  to  pro- 
duce two  vertically  displaced  antenna  beams,  the  vertically  polarized 
wave  for  a similar  configuration  would  result  in  beams  displaced  in 
azimuth.  Since  the  significant  angle  of  arrival  fluctuations  are 
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in  the  vertical  plane,  there  remains  the  question  of  how  to  produce 
vertical  displacement  with  vertically  polarized  waves.  Study  of 
other  modes  used  in  single  horn  monopulse  radars,  e.g.,  and 

TM^^,  has  been  completed  with  only  negative  results.  The  inability 
of  this  technique  to  produce  vertically  displaced  beams  for  both 
polarization  is  a major  limitation. 

A dielectric  lens  can  be  used  either  in  an  end-fire  feed- 
horn  configuration  [4.5]  or  as  a wave  bending  surface  to  produce 
beams  closer  together  for  a fixed  feedhorn  separation.  The  major 
difficulty  of  this  approach  is  the  operation  with  large  power  at 
the  lower  frequencies  required  for  troposcatter  transmission. 

The  use  of  physically  close  feedhorns  to  produce  beams  over- 
lapping at  approximately  the  half  power  beamwidth  is  the  most  attrac- 
tive of  these  approaches  for  reducing  beam  separation.  A further  conr 
sideration  in  the  electronic  steering  of  the  beam  given  adequately  clos 
beams  is  the  requirement  for  beam  steering  of  the  transmit  power. 

Some  of  the  alternatives  include  the  use  of  either  ferrite  phase 
shifters  or  shuttle  switches  to  produce  amplitude  and  phase  differences 
between  beams  to  "point"  the  antenna.  Ferrite  phase  shifters  are 
subject  to  arcing  at  10  kW  power  levels  and  would  require  an  ex- 
tensive waveguide  pressurization  system.  The  shuttle  switch  which 
diverts  a fraction  f of  the  power  to  one  feedhorn  and  1-f  to  the 
other  feedhorn  is  limited  in  state  of  the  art  applications  to  about 
2 kW.  We  conclude  that  although  it  may  be  possible  to  place  the 
beams  within  one  beamwidth  of  each  other  (a  desirable  quality  for 
angle  diversity  also) , electronic  steering  of  up  to  10  Kw  of  power 
appears  to  have  higher  risk  and  cost  than  mechanical  steering 
by  feedhorn  movement. 


4.3  Angle  Diversity 


In  the  last  subsection  we  established  that  vertically  stacked 
feedhorns  located  at  the  focus  of  a parabolic  antenna  produce  an- 
tenna beams  whose  centerlines  are  separated  by  at  least  one  beam- 
width.  For  forward  scatter  systems  one  would  intuitively  expect 
the  energy  received  at  each  feedhorn  to  be  largely  decorrelated 
with  energy  from  other  feedhorns  since  the  respective  common 
volumes  are  almost  disjoint.  Experimental  and  theoretical  studies 
have  confirmed  thst  the  correlation  between  vertical  angle  of  ar- 
rival signals  is  small  enough  to  realize  significant  diversity 
advantage.  If  the  correlation  coefficient  is  below  0.6,  most  of 
the  diversity  is  realized.  The  combining  of  received  angle  of 
arrival  signals  as  separate  diversity  signals  is  termed  angle 
diversity  (AD).  In  forward  scatter  systems,  AD  has  the  distinct 
advantage  over  steering  techniques  in  that  the  diversity  gain  in 
dB  usually  far  exceeds  the  improvement  realized  by  optimally 
pointing  the  antenna.  In  an  AD  system,  the  short  term  fading 
effects  are  mitigated  by  the  diversity  combining  action. 

The  correlation  coefficient  between  beams  is  inversely 
related  to  the  width  of  the  angular  power  spectrum  as  seen  by 
the  antenna  because  of  the  Fourier  transform  relationship  between 
the  spatial  correlation  function  of  the  received  signal  at  the 
antenna  and  its  angular  spectrum.  For  narrow  beam  antennas,  the 
angular  spectrum  in  the  azimuthal  direction  is  limited  by  the 
beamwidth  of  both  antennas  whereas  the  angular  spectrum  in  the 
vertical  direction  is  limited  by  the  beamwidth  of  the  receiving 
antenna  and  the  scattering  angle  dependence.  Since  the  scattering 
angle  dependence  falls  off  more  slovly  than  the  beamwidth  in  most 
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practical  applications,  the  vertical  angular  spectrum  is  broader 
than  the  azimuthal  angular  spectrum.  Thus,  the  correlation 
coefficient  between  vertical  beams  is  smaller  than  for  azimuthal 
beams  given  the  same  displacement.  The  smaller  correlation  coef- 
ficient for  vertical  angle  diversity  beams  makes  it  a superior 
choice  over  horizontal  angle  diversity. 


The  angle  of  arrival  diversity  can  be  realized  by  either  multiple 
transmit  beams,  receive  beams,  or  both.  Multiple  transmit  beams  have 
the  distinct  disadvantage  of  requiring  either  additional  costly  power 
amplifiers  or  a splitting  of  power  into  multiple  feedhorns.  The  3 
dB  loss  in  median  signal  due  to  transmitter  power  splitting  does 
not  result  if  only  one  transmit  beam  aimed  approximately  at  the 
horizon  is  used  and  multiple  receive  beams  are  used  to  intercept  it. 
The  low  correlation  of  vertically  displaced  beams  and  the  power  ad- 
vantage of  multiple  receiver  beams  make  receiver/vertical  angle 
diversity  the  most  attractive  angle  of  arrival  diversity  configuration. 

Long  term  variations  in  median  signal  strength  as  a function 
of  angle  of  arrival  are  compensated  for  in  a steering  system.  How- 
ever, an  AD  system  also  provides  protection  against  long  term  dif- 
fraction conditions  because  of  the  presence  of  additional  angle 
of  arrival  signals.  This  protection  results  in  a long  term  decor- 
relation advantage  of  AD  over  space  or  frequency  diversity  system. 
Because  these  diversity  systems  use  the  same  common  volume,  they 
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are  subject  to  long  term  fading  resulting  from  conditions  effecting 
that  common  volume.  Angle  diversity,  on  the  other  hand,  employs 
virtually  separate  common  volumes  which  in  the  long  term  have 
decorrelated  statistics.  For  example,  earth  bulge  effects  which 
may  seriously  degrade  space  or  frequency  diversity  systems  would 
have  much  less  effect  on  an  angle  diversity  system.  A rriechani- 
cally  steered  system  in  a troposcatter  system  may  provide  a small 
dB  advantage  by  accurately  tracking  the  long  term  optimum  pointing 
angle.  An  AD  system,  because  of  the  beam  separation,  may  not  have 
the  resolution  to  track  small  long  term  variations  but  protects 
against  degradation  due  to  larger  changes  in  angle  of  arrival. 

Thus  in  a forward  scatter  system,  AD  provides  protection  agains 
short  term  fading  effects  and  compensates  for  long  term  gross  varia- 
tions in  angle  of  arrival.  Note  that  in  a transhorizon  diffraction 
radio  system  since  there  is  no  scattering  and  thus  no  multipath 
fading,  it  is  necessary  only  to  track  the  optimum  pointing  angle 
(no  diversity  is  necessary) . If  the  antenna  beam  centerlines  can 
be  separated  by  1 beamwidth,  the  angle  diversity  system  on  a dif- 
fraction path  becomes  a highly  effective  electronic  steered  system. 
As  indicated  in  Section  4,2  however,  practical  alternatives  to  ac- 
complish this  end  are  limited  due  to  the  difficulty  of  electronic 
steering  of  large  transmit  powers. 

4,3.1  Angle  Diversity  Applications 

Because  of  the  diversity  advantage  of  an  AD  system  on  a for- 
ward scatter  link,  the  remainder  of  our  comments  on  angle  diversity 
will  be  directed  toward  forward  scatter  operation.  The  geometry 
for  a vertical  angle  diversity  system  is  shown  in  Fig,  4,5.  Define 
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Fig.  4.5  Angle  Diversity  Geometry 
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Because  of  the  increased  scattering  angle  associated  with 
elevated  beams,  the  received  pol^7er  for  these  beams  will  be  reduced 
and  the  multipath  delay  spread  will  be  increased.  The  effect  of 
the  reduced  power  is  to  limit  the  usefulness  of  the  angle  diversity 
branch  to  protection  against  fades  deeper  than  the  relative  power 
for  that  branch.  The  increased  multipath  will  produce  more  inter- 
modulation distortion  in  analog  systems  and  will  require  equaliza- 
tion in  high  speed  digital  systems.  Horizontally  displaced  feed- 
horns  also  increase  the  scattering  angle  resulting  in  increased 
loss.  The  effect  is  more  promounced  on  longer  paths  than  on 
shorter  paths.  For  shorter  paths  however,  the  correlation  coeffi- 
cient is  large  enough  to  limit  the  diversity  capability. 


In  order  to  assess  the  links  for  potential  application  of 
angle  diversity  we  consider  a simplified  performance  analysis.  The 

improvement  due  to  angle  diversity  can  be  estimated*  by  considering 


These  simplified  calculations  are  presented  here  to  briefly 
outline  the  general  capability  of  angle  diversity  systems. 

In  Section  3,  the  analysis  is  presented  for  the  more  detailed 
and  exact  performance  computations  which  will  follow  in  a 
later  phase  of  the  contract. 


1 


the  average  error  probability  for  DPSK  on  a fading  channel  with  Mth 
order  diversity,  viz., 

- 1 . \-M 


P = - (l  + e^/nJ- 


(4.9) 


The  quantity  is  the  average  received  energy  per  bit 

divided  by  the  noise  power  in  1 Hz.  Equation  (4.9)  is  a good  ap- 
proximation to  performance  of  practical  systems  when  the  implicit 
diversity  effect  due  to  multipath  is  not  large.  If  an  angle  diver- 
sity system  had  D feedhorns  with  independent  relative  received  powers, 
1,  r^,  r^,...r^,  the  average  error  probability  is 

1 S / - 

'’d  Wo)  • 

The  relative  powers  of  the  angle  diversity  branches  must  be 
computed  using  an  angle  diversity  model  whose  theoretical  develop- 
ment is  detailed  in  Sections  2 and  3.  In  order  to  parametrically 
evaluate  the  effects  of  reduced  signal  power  on  the  angle  diversity 
branches,  we  conservatively  estimate  the  relative  powers  by  using 
an  exponential  dependence  on  scattering  angle.  This  estimate  has 
been  shown  in  experimental  studies  to  be  conservative  because  it 
ignores  the  contribution  of  both  the  lower  sidelobe  and  the  lower 
half  beam  width  of  the  elevated  beam. 

The  scattering  angle  dependence  is  taken  as 

1(B)  = sin^^^^e  (4.11) 

For  the  first  elevated  beam,  the  dB  loss  is  conservatively 

estimated  as  1 1 /-s 

sin'^'^O  +a)  P +a 

~ I 11/3  ~ ~ 36.66  log  ^ (4.12) 

sxn  (6_)  o 

o 
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.We^have  determined  the  minimum  practical  beam  separation  to 
be.  one  beamwidth.  Thus 


= 36.66  log(l  + Q/e) 


(4.13) 


which  is  plotted  in  Fig.  4.6. 
scattering  angle  ratio  of 
.o 


0.58 

1.82‘ 


= 0.319. 


The  RADC  link  has  a beam  width  to 


which  gives  a relative  loss  of  4.4  dB. 

In  detailed  AD  calculations  for  the  RADC  path,  we  have 
computed  a smaller  relative  paver  loss  but  the  received  signals 
from  the  two  feedhorns  are  not  independent  as  assumed  in  the  simple 
formulation  presented  here.  The  rather  conservative  estimate  of 
relative  loss  suggested  by  (4.11)  is  balanced  somewhat  by  the 
assumption  of  zero  correlation  between  received  signals.  Thus  the 
simple  parametric  calculation  to  determine  the  general  parameters 
of  troposcatter  circuits  amenable  to  AD  techniques  is  probably  quite 
accurate. 

.As  a general  rule  if  the  elevated  beam  has  a relative  power 
greater  than  one-tenth  it  will  be  useful  contributor  to  a quadruple 
diversity  system.  Weaker  relative  powers  are  useful  in  a dual 
diversity  configuration  where  the  fades  are  deeper.  For  a quadruple 

— _4 

diversity  system  the  error  rate  is  reduced  by  (1+0.1  E^/N^)  for 

a one-tenth  relative  power  AD  branch  per  explicit  diversity  branch 
— -5 

If  P = 10  without  AD,  the  improvement  with  AD  from  (4.9)  and 
(4.10)  is  a factor  of  33  better.  At  2.5  dB/decade  this  corresponds 
to  an  improvement  of  approximately  3.8  dB.  From  Fig.  4.6  for  = 

10  dB  one  finds  that  the  10  dB  contour  corresponds  to 
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n/ 6 = 0.875,  10  dB  controur, 


(4.14) 


i.e.,  AD  is  useful  when  0/ 0 is  less  than  0.875. 

The  expressions  for  bearawidth  and  scattering  angle  are 
60X 

n = 57. 3D  (4.15) 

6 = d/a.  (4.16) 

Prom  these  relationships,  the  10  dB  contour  as  a function 
of  frequency  and  path  length  can  be  computed  and  they  are  plotted 
in  Fig.  4.7.  For  AD  systems  with  length/frequency  parameters  above 
the  contour  a significant  performance  improvement  will  result. 

Note  that  most  DCS  forward  scatter  systems  lie  above  these  contours. 

4.3.2  Summary  of  Previous  Adaptive  Antenna  Control  Experiments 

In  the  field  of  tranhorizon  radio  communication,  empirical 
results  from  radio  experiments  have  continued  to  provide  major 
contributions  to  the  growth  of  the  art.  The  approach  to  an  adaptive 
antenna  control  system  should  be  developed  in  light  of  the  wealth 
of  past  angle  of  arrival  radio  experiments.  In  this  subsection  we 
review  some  of  the  more  significant  of  these  experiments  and  sum- 
marize the  major  findings  to  date  concerning  this  communications 
mode . 

In  an  early  study,  Crawford,  Hogg,  and  Kummer  [4.1]  performed 
a comprehensive  experimental  study  of  transhorizon  tropospheric 
propagation.  The  experiments  included  antenna  pointing  studies 
and  an  investigation  of  angle  diversity.  They  established  that, 
on  the  average,  the  maximum  signal  level  is  received  when  the 
antennas  are  aimed  along  the  great  circle  between  transmitter  and 
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receiver  location,  with  each.being  aimed  just  above  the  horizon. 
However,  they  found  that  local  atmospheric  refraction  can  exert 
considerable  influence  on  the  signals  received  on  narrow  beam 
antennas.  In  their  angle  diversity  experiments  they  found  that 
vertically  displaced  feeds  provided  essentially  uncorrelated  chan- 
nels at  both  460  and  4110  MHz.  Their  results  with  horizontally  dis- 
placed feeds  revealed  a fade  rate  dependency.  The  theoretical 
diversity  improvement  was  obtained  in  almost  all  cases  at  4 GHz  ex- 
cept when  the  fade  rate  was  extremely  slow,  indicating  the  presence 
of  very  large  horizontal  layers.  At  460  MHz  the  theoretical  diver- 
sity improvement  was  not  realized  in  many  cases.  VJhen  the  rms 
Doppler  spread  was  below  approximately  0.5  Hz  (which  it  will  fre- 
quently be  at  this  frequency  of  operation) , large  correlation 
between  the  horizontally  displaced  signals  significant3.y  degraded 
this  type  of  diversity.  In  another  study  of  azimuthal  angle  di- 
versity, Chisholm,  Rainville,  Roche,  and  Root  T4.2]  measured  cor- 
relation coefficients  of  0.7  with  0.8  beamwidth  separation  and 
0.5  with  1.2  beamwidth  separation  on  a 2 GHz,  188  mile  link. 
Significant  diversity  advantage  is  realized  even  with  this  high 
degree  of  correlation.  The  correlation  coefficient  dependency  on 
fade  rate  and  resulting  loss  in  diversity  gain  was  not  investigated 
experimentally.  Vogelman,  Ryerson,  and  Bickelhaupt  [4.3]  measured 
negligible  correlation  between  vertically  displaced  beams  with 
separation  of  approximately  3 beamwidths.  This  large  separation 
leads  to  excessive  power  loss  in  the  elevated  beams,  however,  due 
to  an  increase  in  the  scatter  angle.  The  theoretical  relationship 
between  correlation  coefficient  and  angular  separation  in  beamwidth 
has  been  derived  [4. 3, 4.4]  and  is  shown  in  Fig.  4.8^  Since  con- 
siderable correlation  (<  0.6)  will  not  significantly  degrade  the 
realizable  diversity  gain,  it  is  advantageous  to  place  the  beams 
closer  together  to  reduce  the  power  loss  in  elevated  beams. 
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Pig.  4.8  Correlation  Coefficient  vs.  Angular  Separation 
in  Beamwidths 
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Angle  diversity  te.rts  at  8 GEz  have  been  reported  on  by 

Surenian  [4.5].  A significant  feature  of  this  work  was  the  use 
of  dielectric  material  at  the  feedhorn  opening  to  produce  and  end- 
fire  array.  This  approach  provided  decorr elated  vertical  and 
azimuthal  beams  with  only  a small  median  difference  in  signal 
strength  beams.  There  is  question,  however,  whether  this  approach 
would  be  successful  at  lower  troposcatter  frequencies  because  of 

(1)  increased  dielectric  loss  and  associated  problems  of 
high  power  transmission,  and 

(2)  the  effects  of  slow  fading  conditions  v;hich  increase 
the  correlation  coefficient  between  horizontally  dis- 
placed beams. 

An  angle  diversity  system  was  implemented  for  test  purposes 
on  the  DCS  troposcatter  link  between  Greenland  and  Iceland.  This  ■ 
link  is  characterized  by  800  MHz  FDM/PM  operation,  120-foot  parabolic 
reflectors,  and  a path  length  of  440  miles.  The  test  results  have 
been  reported  by  Travis  [4.6].  They  found  considerable  correlation 
and  relatively  weak  signal  strengths  in  the  horizontally  displaced 
beams  but  very  low  correlation  and  relative  median  signal  strengths 
of  -5,  0,  -2,  and  -8  dB  in  the  four  center  vertically  displaced 
horns.  The  tests  included  the  use  of  an  equal  gain  predetection 
combiner  [4.7]  of  the  type  to  be-de scribed  in  Section  5.  Because 
excessive  multipath  conditions  it  was  necessary  to  modify  the  com- 
biner by  measuring  intermodulation  noise  and  switching  out  diversity 
channels  where  the  combined  thermal  plus  intermodulation  noise  was 
6 to  10  dB  higher  than  the  noise  in  the  best  channel. 

These  system  tests  led  to  the  development  by  Bell  Laboratories 
of  an  angle  diversity  test  system  from  Monrovia,  Maryland  to 
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Crawford  Hill,  New  Jersey.  The  results  of  this  program  have  been 
reported  by  Monsen  [4.8],  Bailey  [4.9],  Merrill  [4,10],  and  Brady 
and  Lima  [4.11]. 


i 
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Some  of  the  more  significant  findings  of  the  program  include: 

• a beam  swinging  experiment  established  that  the  antenna 
pointing  angle  should  be  reduced  roughly  1/2  beaiawidth 
for  each  additional  vertical  feedhorn  (simplex  link) . 

• A triple  angle  diversity  system  was  shown  to  have  ap- 
proximately the  same  outage  probability  vs.  received 
signal  level  as  a conventional  dual  diversity  system. 

• Angle  diversity  exhibits  a long  term  correlation  advan- 
tage due  to  the  long  term  independence  of  the  separate 
common  volumes. 

• The  multipath  delay  spread  increased  by  a factor 
of  1.6  for  the  first  elevated  beam  and  by  a factor 
of  2.7  for  the  second  elevated  beam. 

• Double  humped  multipath  profiles  were  measured  fre- 
quently, particularly  on  the  elevated  beanos. 

• The  Merrill  model  for  prediction  of  angle  diversity 
relative  signal  strength  was  very  accurate  for  the 
second  skybeara  but  somewhat  optimistic  for  the  third. 

• AD  performance  improvement  close  to  the  theoretical 
was  shown  in  0.5  and  1.0  Mb/s  PSK  data  tests. 


Adaptive  antenna  control  with  application  to  knife-edge  dif- 
fraction links  has  been  investigated  in  tests  on  a DCS  link  in 
southern  Europe.  Multiple  feedhorns  were  used  to  determine  if 
outage  conditions  in  the  main  beam  occurred  independently  from  out- 
ages in  elevated  beams.  During  the  test  period  no  instances  of 
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stronger  elevated  beams  were  noted  C4.12].  In  later  tests  [4.13] 
order  wire  control  of  both  transmitter  and  receiver  antennas  was 
used  to  optimize  performance.  Typical  daily  variations  in  diffrac- 
tion conditions  resulted  in  improvements  on  the  order  of  1 or  2 dB 
due  to  optimized  antenna  pointing.  During  storm  conditions,  how- 
ever, more  extreme  refractive  index  changes  led  to  larger  improve- 
ments under  adaptive  antenna  control. 
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SECTION  5 


SYSTEM  CONSIDERATIONS 


The  different  adaptive  antenna  control  approaches  must  be 
considered  with  respect  to  the  type  and  character  of  transhorizon 
radio  systems  within  the  Defense  Communications  System.  This 
section  examines  the  impact  of  AAC  applications  in  the  DCS  environ- 
ment and  proposes  significant  areas  to  be  studied  and  resolved 
during  the  contract  phase.  We  also  address  the  effect  of  time 
delay  differences  in  either  steering  or  angle  diversity  systems. 

The  section  concludes  with  a presentation  on  feasible  combiner  struc- 
tures for  both  analog  and  digital  systems. 

5.1  AAC  Application  for  Diffraction  and  Troposcatter  Paths 

Transhorizon  radio  links  within  the  DCS  utilize  either  forward 
scatter  or  knife-edge  diffraction.  Although  there  are  significantly 
more  scatter  circuits  within  the  DCS,  the  use  of  AAC  techniques  is 
to  be  considered  for  both  modes  in  order  not  to  limit  network  availa- 
bility due  to  outages  on  diffraction  systems.  In  the  forward  scatter 
mode,  diversity  must  be  provided  to  compensate  for  the  short  term 
(0.1  to  5 Hz  fade  rates)  multipath  fading.  Outages  on  these  cir- 
cuits are  produced  by  these  short  term  fading  conditions  in  com- 
bination with  depressed  median  signal  conditions.  The  depressed 
median  conditions  are  most  common  in  winter  months  and  are  largely 
due  to  substandard  refractive  gradients  which  increase  the  scattering 
angle  and  thus  the  path  loss.  Outage  conditions  can  be  mitigated, 
therefore,  most  efficiently  by  adding  additional  diversity.  If 
this  additional  diversity  is  in  the  form  of  angle  of  arrival 
diversity,  the  adaptive  combining  used  in  this  approach  results  in 
a broad  optimum  for  the  antenna  pointing  angle  at  the  receiver. 
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At  the  transmitter,  depressed  median  conditions  due  to  substandard 
refractive  index  gradients  are  best  alleviated  by  pointing  the 
antenna  on  or  slightly  below  the  horizon.  Thus  the  requirement  for 
adaptive  antenna  pointing  in  a receive  angle  diversity  troposcatter 
system  does  not  exist.  Some  of  the  theoretical  considerations 
supporting  this  argument  are  detailed  in  Section  6 and  in  Section  7 
empirical  data  from  a Bell  Laboratories  angle  diversity  test  is  sum- 
marized to  establish  the  existence  of  the  broad  optimum  in  antenna 
pointing  angle. 

For  diffraction  links  the  reverse  situation  is  encountered. 
Propagation  via  the  diffraction  mode  results  in  reception  of  a steady 
nonfading  signal  v^hich  at  worst  may  have  a small  fading  component 
which  will  result  in  Rician  fading.  If  the  fading  component  is 
absent,  no  diversity  at  all  is  required  and  if  it  is  present,  switched 
diversity  such  as  that  used  in  land  based  LOS  systems  or  dual  diversity 
should  provide  adequate  protection.  Outage  conditions  on  diffrac- 
tion links  are  produced  primarily  not  by  short  term  fading  effects 
as  in  troposcatter, but  by  long  term  changes  in  the  refractive  index 
gradient  which  produce  antenna  beam  blockage  by  the  earth's  surface. 
This  effect  is  best  compensated  for  by  adjustments  in  the  transmit 
and  receive  antenna  pointing  angles.  A steered  system  as  opposed  to 
an  angle  diversity  approach  is  therefore  appropriate  for  the  diffrac- 
tion links  within  the  DCS. 

Because  of  the  difficulty  of  electronically  steering  large 
transmitter  powers  (~  10  Kw) , mechanical  steering  of  the  antenna 
pointing  angle  is  most  attractive.  RF  Systems,  Inc.,  is  developing 
a design  of  a weather ized  moveable  feedhorn  structure  for  a typical 
DCS  antenna  configuration.  This  approach  is  simple,  low  cost,  and 
can  be  easily  implemented  within  a closed  loop  adaptive  system  for 
optimizing  the  antenna  pointing  angle. 
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For  troposcatter  systems,  the  angle  diversity  approach  will 
be  reduced  to  practice  and  evaluated  on  the  RADC  test  circuit. 
Within  the  study  phase  of  the  contract,  the  troposcatter  systems 
within  the  DCS  will  be  classified  into  groups  having  approximately 
equal  angle  diversity  characteristics.  Performance  predictions 
for  dual  and  quadruple  diversity  operav.ion  with  and  without  addi- 
tional angle  diversity  will  be  developed  to  assess  the  performance 
improvement  realizable  within  the  DCS.  The  RADC  tests  will  provide 
an  important  empirical  comparison  to  validate  the  prediction  ap- 
proach. 

5 . 2 Time  Delay  Effects 

In  both  steering  and  angle  diversity  systems,  the  larger  chan- 
nel time  delay  associated  with  larger  antenna  pointing  angles  must 
be  accommodated.  The  geometry  for  calculating  the  increase  in  mean 
channel  time  delay  for  an  increase  in  pointing  angle  6 is  shown  in 
Fig.  5.1.  The  chord  of  length  X6  forms  an  approximate  90*^  angle 
with  the  line  segment  X such  that 

a = 90-e 

R = 0. 


Thus 


X5  ^ X6 

sine  e 


(5.1) 


x.iT  = z cose  = z (i-eVz).  (5.2) 

The  path  length  differential  is 

Z+Y-X  = Z-Z  (1-0^2)  = Z9V2 

2 

= x60/2  =d66/4=  d 6/4a,  a = effective  earth  radius. 

(5.3) 

For  a 4/3  earth  radius  the  differential  delay  in  nanoseconds 
is 
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Geometry  for  Calculation  of  Differential  Delay 
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Differential  Delay  = ^Ys2liu57.3)  ' ‘^mi^deg. 

Thus  a 1°  change  in  the  antenna  pointing  angle  on  a 200  mile 
path  will  produce  approximately  176  nanoseconds  of  additional  de- 
lay. The  Megabit  Digital  Troposcatter  Subsystem  (MDTS)  is  designed 
to  operate  at  1.5, 3. 1,6. 3, 9.4,  and  12.6  Mb/s.  At  the  highest  data 
rates,  the  symbol  length  is  160  nanoseconds  and  thus  steering 
adjustments  in  the  antenna  pointing  angle  must  be  made  small 
compared  to  the  symbol  length  to  avoid  loss  of  bit  synchroniza- 
tion. Since  the  steering  adjustments  compensate  for  a long 
term  effect,  the  adjustments  can  be  made  with  steps  which  corres- 
pond to  a differential  time  delay  of  a small  fraction  of  a 
symbol  interval.  For  angle  diversity  both  analog  and  digital  sys- 
tems v/ill  require  delay  padding  of  the  shortest  angle  of  arrival 
path  in  order  to  compensate  for  the  differential  delay.  Since  the 
mean  path  delay  for  troposcatter  and  diffraction  systems  does  not 
vary  appreciably  with  time,  an  initial  delay  pad  insertion  for  each 
feedhorn  would  eliminate  differential  delay  effects  from  both  an 
electronically  steered  system  or  an  angle  diversity  system. 

In  angle  diversity  systems,  one  must  consider  the  method  of 
diversity  combining.  Analog  and  digital  systems  will  be  treated 
separately. 

5.3  Diversity  Combining/Analog  Systems 


(5.4) 


The  combining  of  independent  FDM/FM  signals  can  be  accomplished 
by  a well  known  phase  stripping  technique  [5.1].  The  basic  concept 
is  illustrated  in  Fig.  5.2.  Using  complex  notation,  the  ith  diversity 
signal  has  some  RF  phase  angle  After  multiplying  the  diversity 

signal  by  a reference  signal  and  removing  the  modulation  by  averag- 
ing, the  received  phase  angle  can  be  stripped  off  by  another 
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multiplication  operation.  The  early  combiners  [5.2]  used  a local 
oscillator  to  generate  the  reference  signal  and  mixers  for  the  multi 
plication  operations.  For  this  combiner  the  transmitted  signal  must 
contain  a carrier  component  for  proper  operation.  The  use  of  mixers 
instead  of  multipliers  reduces  the  combiner  from  maximal  ratio  to 
equal  gain  for  which  there  is  not  a large  performance  difference. 
More  recent  combiners  use  a gain  controlled  version  of  the  combiner 
output  as  the  reference  signal  to  produce  what  is  termed  a regener- 
ative predetection  combiner.  This  type  of  combiner  is  stable  and 
does  not  require  a carrier  component.  The  principle  of  operation 
can  be  understood  from  solution  of  an  eigenvalue  problem.  For  the 

regenerative  predetection  combiner,  the  output  can  be  considered  the 

j 9 i 

dot  product  of  diversity  multipliers  w^  v;ith  diversity  inputs  xe  ^ 

y = 2 w.  xe  = w'x  (5.5) 

i=l  ^ 

where  xe  is  the  ith  component  of  the  vector  x. 

The  reference  is  Gy  and  from  Fig.  5.2,  the  weight  vector 
w must  satisfy  the  relationship 
Gy*x  = w 


, 1 
XX  w = — w 


(5.6) 


which  is  an  eigenvalue  problem.  The  regenerative  predetection 

combiner  operates  in  an  identical  manner  to  the  recursion  method 

[5.3]  for  finding  the  largest  eigenvalue  and  its  eigenvector  for 

the  Hermetian  matrix  x easy  to  see  that  the  largest 

eigenvalue  is  1/G  and  its  eigenvector  is  proportional  to  x.  Thus 

in  the  steady  state,  the  diversity  multiplier  weight  vector  is 

proportional  to  x which  both  strips  off  the  RF  phase  and  maximal 

j3i 

ratio  combines  the  signals  i.e.,  each  diversity  input  xe  is 

• * *1  ft  ' 

multiplied  by  the  complex  conjugate  of  the  weight  w^  = xe*^  In 
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the  absence  of  multipath  conditions  the  regenerative  predetection 
combiner  is  optimu.  Although  channel  multipath  in  analog  systems 
introduces  intermodulation  distortion, in  well  designed  analog  sys- 
tems the  bandwidth  is  kept  small  to  avoid  these  distortion  effects. 
This  also  insures  that  the  combiner  will  not  fail  due  to  m.ultipath 
effects.  Under  extreme  multipath  conditions,  intermodulation  dis- 
tortion and  combiner  degradation  will  severely  limit  analog  system 
performance. 

5.4  Diversity  Combining/Pigital  Systems 

In  digital  systems  the  bandwidth  is  usually  large  enouc!-  that 
multipath  effects  are  an  important  consideration  in  receiver  -'■•tign. 
In  fact,  the  multipath  can  be  used  to  improve  performance  by  using 
the  decorrelation  between  m.ultipath  components  to  realize  an  im- 
plicit diversity  gain.  The  MDTS  modem  uses  a decision  feedback 
equalizer  [5.4]  to  eliminate  intersymbol  interference  due  to  multi- 
path  and  to  coherently  recombine  the  multipath  components.  Diversity 
combining  of  channels  is  accomplished  by  use  of  a three  tap  trans- 
versal filter  for  each  diversity  signal.  In  general  any  diversity 
combining  technique  for  high  speed  digital  systems  must  in  some 
manner  cope  with  the  presence  of  significant  multipath.  The  use 
of  a regenerative  predetection  combiner  without  specia.l  a.i:cention 
to  the  multipath  effect  is  not  possible.  The  possible  iryproaches 
for  combining  of  the  angle  diversity  signals  include  the  same  approaches 
for  reception  of  high  speed  digital  data  over  troposcattar  channels. 
The  most  successful  of  these  techniques  are  the  equalizer  concept 
used  in  MDTS  and  a transmitter  time  gating  approach  [5.5J  with 
matched  filter  detection  over  the  received  symbol  length.  The  time 
gating  technique  is  somewhat  simpler  to  realize  because  only  one 
signal  multiplication  per  diversity  branch  is  required.  The  equal- 
izer technique  requires  two  or  more  (MDTS  uses  three)  signal  multi- 
plications per  diversity  branch.  However,  the  time  gating  technique 
requires  a time  gated  transmit  signal  and  its  multipath  protection 
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is  by  necessity  limited  to  less  than  one  symbol  interval.  Further- 
more the  time  gate  at  the  transmitter  increases  the  system  bandwidth 
by  the  reciprocal  of  the  duty  cycle.  The  three  tap  forward  filter 
equalizer  used  in  MDTS  can  handle  multipath  delays  greater  than  one 
symbol  interval  and  this  feature  is  probably  an  overriding  considera 
tion  in  angle  diversity  because  the  elevated  beams  will  have  signif- 
icantly more  multipath  than  the  lower  beams. 

Another  advantage  of  the  equalizer  structure  is  its  striking 
similarity  to  the  regenerative  perdetection  combiner  in  terms  of 
hard'ware  realization.  The  operations  shovm  in  Fig.  5.2  for  the 
regenerative  predetection  combiner  are  identical  to  the  operations 
required  for  weight  generation  and  signal  multiplication  in  the 
MDTS  equalizer.  The  difference  between  the  two  systems  is  the 
multiplicity  of  taps  for  the  equalizer  and  the  type  of  reference 
signal  used.  For  analog  combining  the  reference  is  a gain  con- 
trolled version  of  the  combiner  output.  For  digital  combining  the 
reference  is  a 4PSK  error  signal  which  is  an  instantaneous  measure 

of  the  difference  between  the  sampled  quadrature  components  be- 
fore decision  making  and  the  resulting  decision.  The  generation 
of  the  error  signal  can  be  accomplished  for  any  linear  modulation 
modem  (FAM, 2PSK,QPSK,  Staggered  QPSK,  MSK)  by  delivering  at  the 
modem  output  at  least  one  quadrature  sampled  signal  from  which  a 
comparator  makes  a decision.  From  this  signal  the  error  signal 
is  formed  by  taking  the  sign  of  the  difference  of  the  input  and 
output  of  a comparator  with  the  same  reference  level  as  used  in 
the  modem.  In  many  modems  the  error  signal  is  already  available 
in  the  form  of  a quality  measure,  e.g.,  pseudo-error  techniques 
for  BER  estimation. 

A combiner  which  could  be  used  for  both  analog  and  digital 
system  combining  is  shown  in  Fig.  5.3.  For  the  analog  system  the 
reference  is  the  combiner  output  and  only  one  tap  per  diversity  is 
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used.  In  the  steady  state,  the  weights  + jv^, 

are  the  optimum  weights  for  a maximal  ratio  combiner.  For 
the  digital  system  the  4PSK  error  signal  is  selected  as  the  ref- 
erence signal  and  multiple  taps,  j = 1,2,...J,  in  a transversal 
filter  configuration  are  used.  In  a typical  application  the  taps 
would  be  spaced  at  one-half  the  symbol  interval  and  the  number  of 
taps  required  would  span  the  multipath  spread.  For  most  DCS  applica 
tions,  3 taps  would  probably  be  sufficient. 

In  the  digital  mode,  the  use  of  an  error  signal  derived  from 
the  modem  leads  to  optimum  combining  of  the  diversity  signals  be- 
cause the  sample  and  hold  operations  in  the  modem  are  included  in 
the  combiner  optimization.  In  almost  any  transhorizon  radio  modem 
this  error  signal  is  either  available  in  equivalent  form  or  can  be 
easily  extracted.  In  the  absence  of  this  error  signal,  the  angle 
and  space  diversity  channels  can  be  combined  in  pairs  to  reduce  8 
digital  diversity  channels  to  four.  The  reference  for  this  type 
of  combining  is  the  difference  output  of  a hybrid.  The  sum  output 
is  the  combiner  output.  This  type  of  combining  is  based  on  the 
power  inversion  principle  used  to  eliminate  strong  jammers.  In 
a power  inversion  system  the  adaptation  is  used  to  cancel  signals 
which  are  identical  except  for  a linear  filtering  operation.  A 
block  diagram  of  this  type  combiner  is  shown  in  Fig.  5.4,  The 
equalizers  adapt  to  drive  the  reference  signal  to  a minimum  which 
means  the  sum  output  of  the  hybrid  is  a combined  signal.  A major 
disadvantage  of  this  type  combiner  is  that  it  attempts  to  align 
the  phase  at  every  point  in  time  rather  than  only  at  the  sampling 
time.  Aligning  the  phase  at  every  point  in  time  is  much  more 
difficult  and  unless  a large  nuiriber  of  taps  are  used,  the  com- 
bining may  be  poor. 
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Fig.  5.4  Power  Inversion  Diversity  Combining 


In  summary,  an  equalizer  combiner  can  be  used  to  conbine  DCS 
FDM/FM  signals  by  operation  in  an  analog  mode  and  it  can  be  used 
to  combine  high  speed  serial  data  in  an  optimum  manner  using  a 
modem  supplied  error  signal  or  if  this  signal  cannot  be  obtained, 
suboptimum  combining  by  power  inversion  can  be  used. 

5 . 5 Description  of  Paths 

The  two  types  of  communication  link  under  consideration  are 
sketched  in  Figures  5.5  and  5.6. 

On  the  diffraction  link,  a physical  obstruction  (mountain 
ridge)  is  within  the  radio  horizon  of  both  the  transmitter  and 
receiver;  the  dominant  propagation  mode  is  diffraction  from  the 
highest  points  of  this  obstruction. 

The  optimum  configuration  of  antennas  for  this  link  involves 
aiming  the  antennas  so  their  directions  of  peak  gain  are  lined 
up  with  the  rays  that  terminate  at  the  diffraction  point.  It 
will  be  noted  that  we  have  sketched  these  rays  as  curved  trajec- 
tories. This,  in  fact,  is  a correct  representation  of  the 
normal  downward-bending  of  the  radio  path  due  to  decrease  in 
tropospheric  refractive  index  as  altitude  increases. 

The  principal  source  of  suboptimum  performance  can  readily 
be  visualized  from  the  figure:  when  the  refractive  index  gradient 

changes  as  a result  of  diurnal,  seasonal  or  weather  frontal 
conditions,  the  curvature  of  the  ray  from  either  antenna  to 
the  diffraction  point  will  change,  so  that  the  takeoff  angle  of 
the  ray  relative  to  the  local  horizon  will  change.  If  the 
antenna  boresighting  is  not  changed  in  a corresponding  manner, 
then  the  ray  will  leave  the  antenna  at  a direction  other  than 
that  of  maximtam  gain  causing  a decrease  in  received  power. 

The  tropospheric  scatter  path  sketched  in  Figure  5.6  is 
significantly  different;  the  propagation  mechanism  consists  of 
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scattering  from  tropospheric  inhomogeneities  that  are  simulta- 
neously visible  by  the  transmitting  and  receiving  antennas. 

Since  the  scattering  mechanism  is  strongly  angle-dependent, 
most  of  the  scattered  energy  will  arise  from  points  that  are  at 
the  minimum  height  permitting  mutual  visibility.  The  optimum 
aiming  point  for  a single  beam  is  correspondingly  as  close  to 
the  horizon  as  possible.  Although  refractive  index  variations 
can  cause  changes  in  performance  of  scatter  circuits,  there  is 
no  steering  technique  for  single  beam  antennas  that  leads  to  any 
significant  improvement  over  horizontal  aiming.  For  example, 
if  substandard  refractive  gradients  exist,  the  propagating  rays 
will  more  closely  approximate  straight  lines  so  that  the  lowest 
scattering  point  is  at  a higher  elevation,  the  scattering  angle 
is  greater, and  the  average  path  loss  increases.  However,  the 
strongest  signal  component  still  arises  from  rays  emanating  at 
a horizon  angle  with  zero  tilt  from  the  antennas. 

With  superrefractive  or  ducting  conditions,  on  the  other 
hand,  there  is  a theoretical  possibility  that  a horizontal  ray 
will  be  bent  downward  so  strongly,  that  it  will  intersect  the 
ground,  and  that  the  rays  connecting  transmitter  and  receiver 
will  originate  from  slightly  elevated  angles  of  emission.  However, 
even  though  these  are  not  at  the  peak  gain  of  the  antenna  pattern, 
the  loss  is  more  than  compensated  by  the  major  increase  in  field 
strength  accompanying  the  change  from  scatter  to  refractive  propa- 
gation . 

The  availability  of  dual  beams  from  the  receiving  antenna 
does  permit  a performance  improvement  for  tropo  paths  even  in 
the  absence  of  refractive  anomalies  through  the  use  of  diversity 
reception  techniques.  Figure  5.7  illustrates  the  general  principle 
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Fiaure  5.7  Geometry  for  Angle  Diversity 
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involved.  The  second  receiving  antenna  beam  is  elevated  from 
the  primary  beam  and  also  intersects  the  transmitted  beam;  the 
scatterers  in  the  additional  common  volume  produce  a random 
collection  of  phasor  contributions  to  the  second  beam  output, 
yielding  an  approximately  independent  fading  amplitude  on  the 
second  beam.  In  principle,  this  permits  effective  doubling  of 
the  order  of  diversity  when  the  appropriate  diversity  combining 
equipment  is  available.  The  four  sources  of  degradation  of  this 
desirable  extra  diversity  improvement  are: 

(1)  Lower  power  level  in  upper  beam  due  to  greater 
scattering  angle, 

(2)  Partial  correlation  of  fading  due  to  beam  overlap, 

(3)  Possible  degradation  of  lower  beam  by  inclusion  of 
feed  for  upper  beam,  and 

(4)  Increased  multipath  dispersion  of  upper  beam. 

There  are  several  specific  problems  that  warrant  analytical 
solution  for  the  diffraction  and  scatter  paths. 

For  the  diffraction  path  the  principal  problem  of  interest 
is  the  determination  of  the  path  elevation  angle  as  a function 
of  refractive  index  variations.  This  can  be  converted  into  a 
description  of  incremental  path  loss  due  to  effective  antenna 
aiming  error  for  the  conventional  unsteered  antenna,  which  in 
turn  measures  the  potential  gain  from  steered  beam  techniques. 

For  the  tropopath  the  analytic  problems  involved  are  those 
of  computing  power  level  and  correlation  of  the  signals  on  the 
two  beams,  and  determining  the  multipath  dispersion  of  each  and 
the  correlation  between  individual  multipath  components.  These 
and  other  system  performance  factors  will  be  considered  in  the 
next  section. 
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SECTION  6 
SYSTEM  PERFORMANCE 

The  proposed  adaptive  antenna  control  system  tracks  the  opti- 
mum pointing  angle  for  diffraction  link  application  and  combines 
additional  angle  of  arrival  diversity  signals  in  a space/angle 
or  frequency/angle  troposcatter  application.  The  performance  of 
this  system  is  of  course  dependent  on  the  link  propagation  mode, 
link  geometry,  and  atmospheric  conditions.  Prediction  of  system 
performance  as  of  function  of  link  and  atmospheric  parameters 
will  utilize  the  propagation  model  developed  in  Sections  2 and  3 and  a 
communication  system  model  which  defines  the  performance  of  the 
FDM/FM  terminal  for  analog  systems  and  high  speed  modem  for  digital 
systems.  In  the  Bell  Laboratories  angle  diversity  test  program 
[6.1],  the  prediction  model  was  shown  to  agree  well  with  the  ex- 
perimental results.  The  analytic  base  for  prediction  calculations 
has  been  extended  from  the  Bell  Laboratories  work  to  include  a 
performance  evaluation  of  high  speed  data  transmission  over  trcpo- 
scatter  channels.  For  this  communications  mode  the  performance 
calculation  requires  an  assessment  of  both  the  angle  diversity  im- 
provement and  the  implicit  diversity  improvement  resulting  from 
coherent  combining  of  decorrelated  multipath  components. 

6. 1 Performance  Improvement  from  Beam  Steering 

Steering  of  the  antenna  beam  to  track  the  optimum  pointing 
angle  will  provide  performance  improvement  over  fixed  beam  antenna 
for  transhorizon  radio  systems  which  use  diffraction  as  the  propaga- 
tion mode.  Experimental  evidence  is  also  presented  here  to  show 
that  a broad  optimum  for  the  antenna  pointing  angle  exists  on  angle 
diversity  troposcatter  systems.  Because  of  this  broad  optimum, 
tracking  of  the  optimum  is  not  deemed  necessary  for  troposcatter. 
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For  both  types  of  communication  paths, we  are  faced  with 
analyses  of  curved  paths  over  a curved  earth,  as  indicated  in 
Figure  6.1  (a).  If  we  perform  a bilinear  fractional  transforma- 
tion of  the  plane  including  the  transmitter,  receiver  and  earth 
center,  we  can  map  the  earth  center  into  the  point  at  infinity, 
the  curved  earth  surface  into  a horizontal  line,  and  the  rays 
into  upward-bending  circles.  This  conformal  map  preserves  the 
scattering  angle  at  ray  intersections,  and  preserves  the  angle 
of  a ray  with  respect  to  the  local  horizontal,  but  simplifies 
the  calculation  of  all  of  these  angles  because  of  the  fixed 
Cartesian  coordinates.  Furthermore,  in  the  nev;  geometry  the 
vertical  coordinate  gives  local  altitude  directly. 

The  second  analytical  simplification  is  the  highly  accurate 
one  of  treating  ray  segments  as  parabolic  segments  rather  than 
circular  arcs.  We  therefore  take: 


X = horizontal  displacement  along  path 
relative  to  midpoint, 

h = altitude  (same  units  as  x) , 

and  describe  rays  by  a formula  of  the  general  type 

2 

h = A (x-x  ) + B (x-x  ) . 
o o 


In  this  expression  x^  is  the  ground  terminal  of  the  ray  segment 
and  A is  the  tangent  of  the  angle  the  ray  segment  makes  with  the 
ground;  for  all  practical  purposes  we  can  take  A to  be  the  angle 
itself,  measured  in  radians. 


The  parabolic  coefficient  B is  given  by 

B = ^ + CVN 


1 
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(a)  ^CTU/\L  G^OM^TRY 


CURVED  RAY  ArTtR 
CONEORMAL  MAPPING 


EARTH  SURFACE  AFTER 
CONFORMAL  MAPPING 


(b)  ‘ CONFORMALLY  - MAPPED  GEOMETRY 


Fig.  6,1.  Conformal  Geometry 


where  R is  the  actual  earth's  radius,  C is  a constant  and  'TN 
is  the  gradient  of  the  refractive  index  which  is  negative  in  the 
standard  and  superrefr active  or  ducting  atmospheres,  and  posi- 
tive in  the  strongly  subrefractive  atmosphere,  it  is  a conveni- 
ence for  analysis  in  arbitrary  units  to  rewrite  B as 


B = 


2 R 


eff 


where  R __  is  an  effective  earth's  radius  that  includes  the 
err 


refraction  effects.  It  will  be  noted  that  in  the  ducting  atmos- 


phere we  will  have  negative  values  of  R 


eff* 


We  can  model  a typical  diffraction  path  as  consisting  of  an 
obstacle  of  elevation  H at  x = 0,  with  transmitter  and  receiver 
at  -D/2  and  D/2  respectively,  where  D is  the  total  path  length. 
Since  the  local  elevation  of  a horizontally-aimed  ray  is 


h = B(x  + D/2) 


the  elevation  of  the  ray  at  path  midpoint  is 


h = B D /4. 
o 


Clearly,  for  the  diffraction  effect  to  result  rather  than  tropo- 
scatter,  it  must  hold  that 


H > B D /4. 

Using  the  4/3  earth  as  typical  for  the  standard  atmosphere,  we  have 
R = 3444  nmi. 


R = 4592  nmi 
eff 


so  that 


B = 1.088  X 10  ^(nmi)  ^ 


and  if 


D = 100  nmi 
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then 


H > 0.2722  nmi  = 504  meters  = 1654  feet. 

For  an  arbitrary  value  of  H and  D,  the  ray  to  the  obstacle  from 
the  transmitter  is  given  by  the  equation 

2 

h = A(x  + D/a)  + B(x  + D/2)  , 
so  that  at  X = 0 we  have 

H = A D/2  + B D^/4, 

and  consequently  the  ray  takeoff  angle  A is  given  by 
A = (2/D)  (H  - BD^4)  . 

The  condition  for  ducting  is  that  the  refractive  gradient  be 
AN  « -157  N-units/km. 

The  normal  refractive  index  gradient  is  about  1/4  to  1/3  of  the 
ducting  value,  corresponding  to  effective  earth  radius  equal  to 
4/3  to  3/2  of  actual  earth  radius.  In  most  locations  on  the  earth, 
gradients  ranging  from 
AN  = +50 
to 

AN  = -150 

are  the  only  values  that  will  be  seen  for  any  appreciable  percent- 
age of  the  time.  The  corresponding  effective  earth  radius  will 
range  from  about  3/4  actual  radius  up  to  approximately  infinity 
while  the  corresponding  B values  will  range  from 

B = 0 (superrefractive) 

-4  -1 

through  B = 1.1  x 10  (nmi)  standard 


to 


-4  -1 

B = 1.9  X 10  (nmi)  , subrefractive. 


From  the  previous  formula  for  takeoff  angle  we  see  that 


min 


max 


2H/D  - 10  ^ D when  B 2 x 10  ^ 


2H/D 


when  B 
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so  that 


-4 

A - A . « 10  D 

max  min 

no  matter  what  the  obstacle  height.  Thus,  the  variation  in  take- 
off angle  for  a knife-edge  diffraction  link  is  porportional  to 
link  distance  and  independent  of  obstacle  height. 

A path  length  of  175  nautical  miles  could  therefore  have 
a spread  of  elevation  angles  of  1°  or  upward  and  down\^ard  fluctua- 
tions of  . Correspondingly  fixed-steered  antennas  with  1°  beam 
width  at  each  terminal  could  experience  a 3 dB  loss  at  each  rer- 
minal,  while  beam  width  antennas  could  experience  a complete 
loss  of  signal.  Thus  for  diffraction  systems,  steering  can  eliminate 
complete  outage  conditions  due  to  changes  in  the  refractive  index 
gradient. 

The  potential  performance  improvement  from  optimum  antenna 
pointing  under  varying  refractive  index  conditions  is  small  for  the 
forward  scatter  mode.  Considerable  empirical  data  to  support  this 
result  has  been  collected  and  evaluated  in  a beam  steering  experi- 
ment conducted  during  the  Bell  Laboratories  angle  diversity  tests[6.l]. 

In  these  tests  three  vertical  feedhorns  were  used  to  receive 
angle  of  arrival  signals  and  antenna  pointing  at  the  receiver  was 
accomplished  by  movement  under  manual  control  of  the  60'  parabolic 
reflector.  The  path  length  was  180  miles  and  a frequency  of  2.1 
GHz  was  used  in  the  tests. 

The  beam  swinging  experiment  consisted  of  33  days  of  data  col- 
lection spaced  between  September  24  and  December  30,  1970.  The 
data  consisted  of  measuring  the  median  received  powers  for  each 
feedhorn  as  a function  of  antenna  pointing  angle.  The  antenna  at 
the  Crawford  Hill  receiver  has  a beamwidth  of  0.6°  and  was  varied 
in  0.2°  steps  between  elevation  angles  of  -0.6°  to  0.8°.  The 
elevation  angle  is  defined  as  the  angle  between  the  beam  center 
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■ ~>x. 


mi- 


of  the  main  fo^edhorn  and  the  horizontal.  The  horizon  ray  was  de- 

o . o 

pressed  by  0.1  so  bores ight  on  the  horizon  corresponded  to  -0.1  . 

The  data  was  reduced  by  calculating  the  outage  probability  for 
dual  and  triple  angle  diversity  for  each  of  the  33  experiments. 

For  each  of  the  experiments  the  optimum  pointing  angle  was  deter- 
mined at  each  outage  level.  Deviation  from  the  optimum  by  as  much 
as  a 1/2  beamwidth  in  either  direction  produced  little  change  in 
the  outage  level.  The  following  table*  v;hich  is  a computer  print- 
out of  the  reduced  data  for  the  October  20  experiment  is  typical 
of  all  the  experiments  in  revealing  the  broad  optimum  for  antenna 
pointing  angle.  Higher  order  diversity  systems  would  sharpen  the 
optimum  by  the  Dth  power  of  the  outage  probabilities  determined. 
The  optimum  still  remains  broad,  however,  and  one  concludes  that 
there  is  little  performance  improvement  to  be  realized  by  adaptive 
steering  of  a vertical  angle  diversity  system. 

6.2  Performance  Improvement  for  Angle  Diversity 


The  performance  improvement  for  angle  diversity  systems  depends 


on  the 


(1)  system  and  link  configuration 

(2)  diversity  configuration 

(3)  received  power  per  diversity  channel 

(4)  correlation  between  diversity  signals  and 

(5)  multipath  structure. 


We  have  outlined  system  trade-off  considerations  which  suggest 
a vertical  angle  diversity  augmentation  to  existing  troposcatter 
systems.  The  elevated  angle  diversity  beams  can  be  derived  from 
placing  feedhorns  below  the  main  beam  feedhorn.  Since  it  is  desir- 
able to  place  the  beams  close  together  in  the  sky,  a replacement 
* o 

Subtract  1 from  the  reference  elevation  readings  to  obtain  the 
elevation  angle  defined  in  the  text. 
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of  the  existing  feedhorn  structure  with  a new  structure  with 
smaller  feedhorns  for  each  beam  has  been  suggested.  In  almost 
all  link  applications,  no  more  than  dual  angle  diversity  will  be 
recommended  as  an  effective  cost  and  performance  system  improve- 
ment. 

Diversity  configurations  are  identified  by  the  number  of  dif- 
ferent diversity  types  and  in  almost  all  applications  the  total 
diversity  order  is  equal  to  the  product  of  the  number  of  diversity 
channels  in  each  type.  Thus  a space  (S) ,N^  frequency  (F) , 
angle  (A) , and  polarization  (P)  system  would  have  a total  of 
^1  ^2  ^3  ^4  channels.  Note  that  the  effective  diversity 

order  may  be  smaller  than  this  product.  Some  of  the  most  typical 
diversity  configurations  are  summarized  in  Table  6.1. 


Diversity 

Type 

2S 

2F 

2A 

2S/2F 
2 S/2  A 
2S/2F/2A 
2S/2P/2A 


Table  6.1 


Diversity  Configurations 
Total 

Diversity  Nuxriber  of 

Order  Power  Amplifiers 


2 1 

2 2 

2 1 

4 2 

4 1 

8 2 

8*  2 


Additional 
Power  Factor 

0 dB 

3 

0 

3 

0 

3 

3 


Note  that  the  system  with  additional  frequencies  or  polariza- 
tion require  an  additional  power  amplifier  to  realize  the  same  re- 
ceived diversity  signal-to-noise  ratio.  Thus,  the  space  and/or 
* 

The  effective  diversity  order  for  this  system  is  approximately  6th 
order  due  to  correlation  of  the  cross  paths  [6.3]. 
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angle  diversity  systems  enjoy  an  initial  3 dB  advantage  for  a 
fixed  transmitted  power.  Since  any  reasonable  performance  com- 
parison should  be  based  on  the  same  transmitted  power  for  the 
compared  systems,  the  relative  dB  advantages  and  performance  curves 
will  be  normalized  to  a 0 dB  power  factor  or  equivalently  a single 
power  amplifier.  One  should  recognize  that  multiple  frequency 
systems  require  two  power  amplifiers  (since  the  amplifiers  are 
operated  in  saturation,  signals  cannot  be  combined  prior  to  final 
amplification)  and  thus  the  3 dB  disadvantage  of  frequency  diversity 
systems  is  reflected  in  additional  cost  rather  than  system  perfor- 
mance margin. 

The  received  pov/er  per  diversity  channel,  the  correlation 
statistics,  and  the  multipath  structure  for  the  performance  analysis 
will  be  determined  using  NBS  prediction  methods  to  obtain  the  main 
beam  path  loss  distribution  and  the  results  of  the  prediction  model 
described  in  Sections  2 and  3 for  computation  of  coupling  loss, 
relative  received  power,  and  the  cross-channel  multipath  profile, 

. This  profile  is  a complete  description  of  the  correlation 
and  statistical  properties  of  the  multi-filter  uncorrelated  scatter 
channel.  Before  discussing  in  greater  detail  the  computations  re- 
quired for  performance  assessment  of  digital  systems,  we  consider 
first  analog  system  performance. 

For  analog  systems,  diversity  improvement  is  measured  by  the 
outage  distribution, i.e. , the  probability  that  a given  fade  depth 
will  be  exceeded.  The  predetection  combiner  outlined  in  Section 
5.3  provides  maximal  ratio  combining  for  which  the  outage  distribu- 
tion can  be  calculated  for  unequal  and  correlated  diversity  branches 
by  a transformation  of  the  complex  Gaussian  received  processes  to  a 
set  of  independent  received  processes.  This  is  a standard  calcula- 
tion for  correlation  diversity  systems  [6.2]  and  the  details  are 
omitted.  The  outage  distribution  for  analog  transmission  over  the 
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C-band  RADC  test  link  configured  for  dual  vertical  angle  diversity 
has  been  computed  as  an  example.  For  a 168  mile  path  with  28  foot 
antennas  at  4.46  GHz,  the  elevated  beam  has  a relative  loss  of 
1.7  dB  and  a power  correlation  coefficient  of  0.52.  Antenna  side- 
lobes  were  not  included  in  this  first  calculation  but  the  computer 
program  is  being  modified  to  incorporate  their  effect.  The  cal- 
culated outage  distribution  for  this  C-band  troposcatter  link  is 
given  in  Fig. 6. 2* 

The  dual  space/dual  frequency  curve  assumes  only  one  power 
amplifier  at  the  transmitter  which  is  all  that  is  required  for  a 
dual  space/dual  angle  system.  The  conventional  quad  diversity 
configuration  uses  two  power  amplifiers  which  also  could  be  used 
in  the  dual/space  dual  angle  configuration  to  improve  performance 
by  3 dB.  The  dual  space/dual  frequency/dual  angle  also  assumes 
one  power  amplifier  and  could  be  improved  by  3 dB  with  the  use  of 
the  second  power  amplifier  present  in  a conventional  configuration. 
The  improvement  of  an  eighth  order  system  using  angle  diversity 
over  a conventional  quad  system  is  seen  to  be  7 dB  for  an  outage 
level  of  0.001%, 


It  is  generally  desirable  to  use  two  power  amplifiers  in  an 
operational  situation  in  order  to  provide  redundancy  in  case  of 
amplifier  failure.  A dual  space/dual  angle  configuration  requires 
the  use  of  two  power  amplifiers  at  the  same  frequency.  One  ap- 
proach to  accomplish  the  end  is  to  add  the  power  amplifier  outputs 
in  a summer  and  direct  the  sum  port  output  to  one  antenna  and  the 
difference  port  output  to  the  other  antenna.  To  reduce  the  scal- 
loping of  the  antenna  pattern  resulting  from  equal  power/displaced 
antennas,  a phase  shifter  at  the  exciter  output  of  one  transmitter 
can  be  used  to  minimize  one  of  the  splitter  outputs. 
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For  digital  systems  the  improvement  due  to  angle  diversity 
must  include  the  capability  of  modern  modems,  e.g.,  MDTS,  DAR-4, 
to  realize  implicit  diversity  from  multipath  contributions.  This 
is  particularly  important  with  angle  diversity  because  the  elevated 
beam  contains  more  multipath.  We  outline  here  a performance  pre- 
diction model  for  the  MDTS  modem  which  will  be  useful  for  predicting 
performance  during  the  proposed  digital  data  tests  on  the  RADC  link 
and  for  projecting  the  angle  diversity  performance  on  DCS  links 
in  general. 

The  performance  prediction  model  is  an  analysis  which  computes 
the  average  BER  as  a function  of  the  ratio  of  average  bit  energy 
to  noise  power  in  1 Hz,  i.e.,  E /N  . The  computation  includes  the 
MDTS  modem  degradations  which  significantly  affect  the  average  BER. 
These  degradations  account  for  approximately  a 2 to  3 dB  shift  from 
the  ideal  CPSK  error  probability  curve.  The  major  contributors 
the  degradations  are  (1)  imperfect  90°  splitters  in  the  transmitter 
and  receiver,  (2)  non  ideal  matched  filtering,  (3)  time  tracker  off- 
set, (4)  and  adaptive  equalizer  misadjustment . The  analysis  con- 
siders the  number  N of  forward  equalizer  taps  per  diversity  chan- 
nel as  the  rank  of  a linear  vector  space.  For  the  complex  Gaussian 
uncorrelated  scattering  troposcatter  model,  an  orthogonal  transforma- 
tion exists  in  this  vector  space  which  converts  the  available  im- 
plicit diversity  and  correlated  angle  diversity  into  independently 
fading  diversity  channels.  Each  of  these  diversity  channels  has 

the  fraction  X.,  i = 1,2,...N  of  the  average  E„/N  associated  with 
1 B o 

it.  The  X^'s  are  the  eigenvalues  for  the  transformation  and  they 
satisfy 

N 

S X . 5 1,  0 < X . , i = 1,2, . . .N. 

i=l  ^ ^ 
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In  a manner  of  speaking  the  tinalysis  converts  the  implicit  diver- 
sity and  correlated  angle  diversity  structure  present  in  the  func- 
tion Q'tm(§)  to  an  angle  diversity  structure  with  each  angle  diver- 
sity path  receiving  independent  power  contributions  i = 1/2,... N 
of  the  average  received  power. 

The  computation  of  the  optimum  equalizer  weights  after  this 
transformation  is  tribial  as  the  equalizer  reduces  to  the  maxi- 
mal ratio  combiner.  Since  the  average  BER  performance  for  the 
latter  is  v;ell  known  for  PSK  systems,  the  orthogonal  transforma- 
tion make  calculation  of  the  average  BER  as  a function  of  implicit 
diversity  possible.  The  orthogonal  transformation  considers  only 
one  bit  and  neglects  intersymbol  interference. 

Because  the  intersymbol  interference  is’ neglected  the  perfor- 
mance calculated  is  a lov/er  bound  which  ignores  the  loss  in  implicit 
diversity  due  to  the  need  to  remove  intersymbol  interference.  This 
loss  which  is  termed  the  intersimribol  interference  penalty  is  less 
than  a dB  for  small  to  moderate  multipath  spreads, e.g. , ratios  of 
twice  the  rms  multipath  spread  {2q)  to  baud  length  (T),  which  satisfy 

2a/T  <0.5. 

For  the  larger  ratios  the  intersyn^ol  interference  penalty  is 
larger  and  the  bound  is  not  tight.  A mathematical  model  to  include 
the  effects  of  intersymbol  interference  and  thus  obtain  performance 
predictions  for  2a/T  > 0.5  has  been  developed.  The  model  validity 
has  been  established  by  comparing  the  model  results  with  MDTS  sim- 
ulator and  empirical  results.  Currently  a computer  program  modifica- 
tion of  the  model  to  include  the  effects  of  angle  diversity  correla- 
tion is  being  undertaken.  The  modification  will  utilize  the  numerical 
integration  outputs  of  the  analytic  procedure  described  in  Sections 
2 and  3 for  determining  the  cross-channel  multioath  profile  Q,  (§)  . 
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I In  the  next  interim  technical  report,  a description  of  the  modem 

^ prediction  technique  for  angle  diversity  channels  will  be  given 

! along  with  performance  results  for  typical  and  actual  troposcatter 

; link  configurations . 
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APPENDIX  A 


SCATTERED  FIELD  SOLUTION 

In  this  appendix  the  solution  from  Maxwell's  equations  of 

the  scattered  field  at  the  receive  antenna  is  obtained  in  terms 

of  the  incident  field  in  the  scattering  volume  and  the  refractive 

index  fluctuations.  The  derivation  is  based  on  the  original  scat- 

* 

tering  field  solution  derived  by  Tatarski  . 

Consider  the  tropospheric  scatter  geometry  shown  in  Fig.  A.l. 
For  an  incident  field  distribution  A(r')  in  the  scattering  volume 
we  wish  to  determine  the  electric  field  at  the  observation  point 
r_  from  the  origin  0 near  the  center  of  the  scattering  volume.  Be- 
cause of  the  distance  between  scattering  volume  and  antennas  we 
can  take  both  the  incident  and  scattered  fields  to  be  plane. 

We  assume  that  the  magnetic  permeability  u is  unity.  The 

2 

permittivity  e and  the  refractive  index  n are  related  by  e = n . 

The  refractive  index  is  composed  of  a fixed  mean  value  n near  unity 
and  a random  component  n^. 

n = n + n^ 

It  is  the  presence  of  the  random  component  which  produces 
scattered  energy  back  to  the  receiver  on  the  earth's  surface.  With- 
out loss  of  generality  we  can  take  the  mean  value  of  the  refractive 
index  to  be  unity  and  account  for  any  deviations  from  this  assump- 
tion by  replacing  the  propagation  constant 


* 


V.I.Tatarski  /'Wave  Propagation  in  a Turbulent  Medium,"  Trans- 
lated by  R. A. Silverman,  McGraw- Kill, New  York,N.Y.,  1961. 


Fig.  A.l  Troposcatter  Geometry 


by  kn  in 
and  thus 

the  final  results. 
Maxwell ' s equations 

We  assume  an  e time  dependence 

in  CGS  units  take  the  form 

7 X E = jkH 

(A.l) 

7 X H = -jkeE 

(A.  2) 

7 • =0  . 

(A. 3) 

The  permittivity  e because  of  its  random  component  must  be 
considered  as  a spatial  variable.  Using  the  identity 

2 

7 X (7  X E)  = - 7 E + 7(7  *E) 
and  the  first  of  Maxwell's  equations,  we  obtain 

2 2 
“7E  + 7(7  *E)  = 7X  (jkH)  = k eE  . (A. 4) 

From  the  divergence  equation 

0=7  • eE=s7  »E+E*  7e  (A. 5) 

or 

v*E  = -—  E*7e  = -E  - 7’tne. 

— e ~ ~ 

Combining  (A. 4)  and  (A. 5)  we  have 

2 2 2 
-7  E - 7(E  *7-tns)  = k e E 

or 

2 2 

7 E + 27 (E  •7vnn)  + k nE  = 0 (A. 6) 

2 

where  we  have  substituted  e = n . We  assume  that  the  fluctuations 
of  the  refractive  index  are  small  compared  to  the  mean,  i.e., 

2 

n^  « 1 

so  that  we  can  make  an  expansion  of  the  electric  field  of  the  form 

E = E + E,  + E„  + ... 

— -TO  —1  —2 
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where 


E.  is  the  order  smallness  of  . 

1 1 

Keeping  only  the  first  two  terms  of  the  electric  field,  and 
using  the  first  order  approximation  -tn(l+n^)=  n^,  the  expression 
in  terms  of  the  first  two  field  component  becomes 

(E  + EJ  + 2V  r(E  + E,  ) -vn,"!  + (l+2n,  ) (E  + Ej  =0. 

— O “1  L ~0  ~1  1 ~0  “1 

If  we  group  terms  of  order  1 and  n^,  we  have 

V^E  + k^E  = 0 (A. 7) 

— o — o 


V^E,  + k^  E,  + 2k^n,E  + 2v(E  • vn)  = 0 
—1  —1  1—0  — o — 


(A. 8) 


The  solution  to  (A. 7)  at  a field  point  r is  the  spherical  wave 


E = A (r)  e 
— o — o — 


jk’r 


which  is  taken  as  the  incident  wave  with  propagation  vector  k=k(r/r) 
The  cartesian  components  of  (A8)  are  of  the  form 

v^u  + k^u  = f (r) 


which  has  the  Green's  function  solution 


1 n 3 

u(r)  = - — J f(r')  dr'  (A. 9) 

^ |r  - r'  I 

v 

where  r'  is  a variable  vector  ranging  over  the  scatter  volume  V. 
Since  the  receive  antenna  is  much  farther  away  than  the  linear 
dimension  of  the  scattering  volume,  the  denominator  in  (A. 9)  can 
be  approximated  by  lrl=r.  The  solution  for  the  scattered  field  is  then 
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% (r) 


Jsl 

2nr 


p _ , , . , , , jk( Ir-r ' I ) 3 

1 E (r ' ) n,  (r ' ) e dr 

j o 1 


V 


+ ^ r ^Iklr-r'l  . 3^,  A ^ 

2nr  J — o 1 12 


V 


The  second  integral  can  be  integrated  by  parts  (using  Gauss's 
Theorem) 


I X vy  = J xy  - I yvx 


V 


V 


and  the  surface  integral  disappears  since  the  volume  can  be  made 
large  enough  that  the  integrand  is  negligible.  Thus 


^2  “ " 2rrr 


(E  • yn,)  d’r' 

~o  1 


ik  p,_  . ik|r  - r'l  ,3  . 

^ — J e ~ ' —T  d r' 


2rrr 


r-r 


V 


= ar  J (E  'Vn  ) e^^^—  — ^ d^r' 

2nr  — J — o 1 

V 

where  a^  is  the  unit  vector  r/r.  Integrating  by  parts  again  and 
neglecting  the  surface  integral  gives 


I = a 
■^2  2TTr  -r 


J ^ 


V - E e 
— o 


jklr-r'I  ^3^, 


v 


and 
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„ „ jklr-r'l  jklr-r'l  ' r,  r,  jkir-r' 

V *E  e ' ' = e ‘ ' V • E + E • ve-^ 


— o 


— o — o 


„ . jk|r-r'|  _ -r+r' 

= 0 + 3ke ' E * 

-o  r-r‘ 


• .-u  Dk  r-r 

= -ike  ‘ 'E  • a 

— o — r 


where  we  have  used  the  fact  that  the  divergence  of  the  incident 
field  is  zero.  The  two  integrals  are  now  seen  to  be  identical 
except  for  their  vector  direction.  The  scattered  rield  is  thus 


E f[E  - a (E  . a )1n,  . 

—1  2^r  JL“o  — r — o — r J 1 

V 

The^ vector  E^  is  easily  seen  to  be  perpendicular  to  the  unit 
vector  Let  a^  be  any  unit  vector  satisfying  ^ = 0,  then 

_ k^  p . jklr-r'l, 3 , 


E.  • a = 0. 
—1  — r 


The  incident  field  E has  a vector  direction  which  makes  an 

— o 

angle  with  the  unit  vector  to  the  source  a of  90°-P . where  P is 

_r 

the  scattering  angle.  For  small  scattering  angles,  the  dot  product 
a^^  * E^  is  approximately  equal  to  the  incident  field  function 

~ . Thus  tne  electric  field  at  the  receiving  antenna 
is  a transverse  plane  wave  with  amplitude 


_ , . f n / f \ ^ i\  jk(|r-r'j  + m*r')  ,3  , 

^1^^^  2^  J "^1  - ® ^ 


V 

which  is  the  expression  used  in  Section  2 to  derive  the  cross- 
channel multipath  profile. 
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